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ABSTRACT 


This  report  covers  two  years  of  work  by  the  Jornada  Experimental  Range  Department  at  New  Mexico 
State  University  and  the  U.S.Army  ERDC/CERL.  The  overall  goal  of  the  project  was  to  enhance  the 
vegetation  simulation  modeling  capabilities  of  Fort  Bliss,  Texas,  for  evaluating  the  effects  of  military 
disturbances  on  arid  land  vegetation.  The  aim  of  this  work  was  to  mechanistically  predict  shifts  from  one 
community  type  to  another  by  determining  primary  linkages  between  community  types. 
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PREFACE 

This  study  was  conducted  for  the  Office  of  the  Directorate  of  Environmental 
Programs  (DAIM),  Assistant  Chief  of  Staff  (Installation  Management) 
(ACS[IM])  under  project  622720896,  “Environmental  Quality  Technology”; 
Work  Unit  CNN-T081.  The  technical  monitor  was  Dr.  Vic  Diersing,  DAIM-ED- 
N. 

The  Construction  Engineering  Research  Laboratory  (CERL)  Principal 
Investigator  was  Dr.  Jeffrey  S.  Fehmi.  The  managers  at  Fort  Bliss  were  Kevin 
vonFinger  and  Brett  Russell.  The  work  was  performed  under  contract  by  Debra 
P.  C.  Peters  and  Tamara  Hochstrasser  in  association  with  New  Mexico  State 
University  and  USDA-ARS  Jornada  Experimental  Range.  The  technical  editor 
was  Gloria  J.  Wienke,  Information  Technology  Laboratory.  Mr.  Stephen  Hodapp 
is  Chief,  CEERD-CN-N,  and  Dr.  John  T.  Bandy  is  Chief,  CEERD-CN.  The 
associated  Technical  Director  was  Dr.  William  D.  Severinghaus,  CEERD-CV-T. 
The  Director  of  CERL  is  Dr.  Alan  W.  Moore. 

CERL  is  an  element  of  the  U.S.  Army  Engineer  Research  and  Development 
Center  (ERDC),  U.S.  Army  Corps  of  Engineers.  The  Commander  and  Executive 
Director  of  ERDC  is  COL  John  Morris  III,  EN  and  the  Director  of  ERDC  is  Dr. 
James  R.  Houston. 


Simulation  of  Vegetation  Recovery  from  Military 
Disturbances  on  Fort  Bliss 


TAMARA  HOCHSTRASSER,  DEBRA  P.C.  PETERS,  AND  JEFFREY  S.  FEHMI 

1  INTRODUCTION 


This  report  covers  two  years  of  work  by  the  Jornada  Experimental  Range 
Department  at  New  Mexico  State  University  and  the  U.S.  Army  ERDC/CERL. 
The  overall  goal  of  the  project  was  to  enhance  the  vegetation  simulation  model¬ 
ing  capabilities  of  Fort  Bliss,  Texas,  for  evaluating  the  effects  of  military  distur¬ 
bances  on  arid  land  vegetation.  The  project,  conducted  by  the  Jornada  Experi¬ 
mental  Range  (JER),  where  research  on  arid  land  vegetation  dynamics  has  been 
conducted  for  90  years,  was  crucial  to  transferring  scientific  knowledge  on 
vegetation  dynamics  to  Fort  Bliss. 

The  aim  of  this  contract  work  was  to  mechanistically  predict  shifts  from  one 
community  type  to  another  by  determining  primary  linkages  between  community 
types.  Both  community  types  and  linkages  between  community  types  can  be 
defined  in  different  ways.  For  the  purposes  of  this  project,  we  use  the  following 
definitions: 

Community  type:  A  community  was  defined  based  on  the  most  important 
plant  species  (>50%  of  total  cover).  A  shift  from  one  community  type  to  another 
was  defined  as  a  shift  in  the  composition  of  plant  species  that  constitute  >50%  of 
the  total  cover  of  the  vegetation.  This  definition  has  become  more  acceptable  in 
ecological  research,  as  the  biotic  structure  of  ecosystems  is  predominantly  due  to 
the  presence  of  plants,  the  majority  of  which  belong  to  the  dominant  plant  species 
(by  definition). 

Linkages  between  community  types:  Community  types  are  linked  by  mortal¬ 
ity  and  recruitment  events  of  dominant  plants.  A  decrease  in  the  population  of 
one  dominant  species  and  an  increase  in  the  population  of  another  one  was 
defined  as  a  shift  in  community  type  as  defined  above  (Fig.  1).  In  arid  environ¬ 
ments,  mortality  and  recruitment  events  of  plants  are  affected  by  water 
availability,  seed  dispersal,  competition  with  existing  vegetation,  and  disturbance. 
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While  the  ultimate  goal  of  this  project  was  to  predict  transitions  between  all 
community  types,  the  primary  focus  was  to  predict  shifts  between  grass-  and 
shrub-dominated  systems.  A  shift  in  the  dominant  lifeform  leads  to  more  signifi¬ 
cant  ecological  changes  in  the  system  than  a  shift  between  dominant  species  of 
the  same  lifeform  (Archer  et  al.  2001).  In  arid  ecosystems,  a  shift  in  the  dominant 
lifeform  affects  water  dynamics  (Schlesinger  et  al.  1 990),  leads  to  a  redistribution 
of  nutrients  in  the  landscape  (Schlesinger  et  al.  1990;  Schlesinger  et  al.  1996), 
and  influences  functional  group  composition  of  subdominant  plant  species 
(Hochstrasser  and  Peters,  submitted).  Furthermore,  the  distribution  and  abun¬ 
dance  of  animals  in  the  ecosystem  may  also  be  affected  (Parmenter  and 
MacMahon  1983;  Whitford  1997;  Pidgeon  et  al.  2001).  Therefore,  disturbances 
that  alter  the  dominance  from  grasses  to  shrubs  (or  vice  versa)  affect  the  system 
more  fundamentally  than  disturbances  under  which  the  system  maintains  the 
dominant  lifeform.  For  the  purposes  of  this  project,  it  was  therefore  important  to 
understand  the  ecological  mechanisms  that  change  the  dominant  lifeform,  and 
how  human-induced  disturbances  can  influence  these  mechanisms. 
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The  three  major  objectives  for  understanding  and  minimizing  the  ecological 
impacts  of  military  training  activities  are  compliance  with  the  law,  maintenance 
of  the  landscape  for  effective  training,  and  reduction  of  off-site  environmental 
pollution  (Goran  et  al.  1983).  The  laws  that  regulate  the  management  of  ecologi¬ 
cal  impacts  on  military  lands  are  the  National  Environmental  Policy  Act  (NEPA), 
Army  Regulations  200-285,  and  the  Endangered  Species  Act.  Military  training  is 
less  effective  if  there  is  a  lack  of  cover  for  people  and  vehicles,  and  if  severe 
environmental  degradation  reduces  the  trainability  of  the  area.  Dust  from  the 
training  area  may  potentially  impact  the  air  quality  of  communities  downwind, 
and  sediment  in  streams  eroded  from  degraded  land  may  affect  water  quality 
downstream. 

Recovery  from  disturbance  in  arid  lands  generally  takes  a  long  time.  In  the 
Mojave  Desert,  recovery  from  soil  compaction  in  abandoned  mining  town  sites 
can  take  over  100  years  (Webb  et  al.  1986).  Mechanisms  for  recovery  of  com¬ 
pacted  soil  could  be  freezing  and  thawing  cycles  and  the  activity  of  termites 
(Taylor  Soltero  1996).  For  managing  arid  lands,  it  is  crucial  to  estimate  recovery 
times  from  disturbances.  Generally,  estimating  recovery  time  of  the  vegetation 
has  been  done  using  linear  extrapolation  (Lathrop  1983;  Webb  et  al.  1986).  How¬ 
ever,  recovery  dynamics  may  not  be  linear  over  time,  and  a  process-based 
approach  to  extrapolation  may  be  more  insightful.  The  goal  of  this  project  was  to 
provide  process-based  approximations  of  recovery  times  from  disturbances  of 
desert  grasslands  and  shrublands. 

A  previous  effort  at  creating  a  simulation  model  for  evaluating  the  effects  of 
military  disturbances  on  vegetation  at  Fort  Bliss  was  the  EDYS  model  (Childress 
et  al.  2002).  This  model  focused  at  the  landscape  scale  but  incoiporated  a  multi¬ 
level  spatial  nesting  structure  from  the  plant  to  the  plot  to  the  landscape.  The 
relative  performance  of  individual  plants  given  a  certain  climate  and  disturbance 
regime  was  evaluated  at  the  plant  scale.  The  composition  of  the  vegetation  at  the 
plot  scale  was  then  translated  into  a  vegetation  type.  Once  the  vegetation  compo¬ 
sition  sufficiently  changed,  the  plot  shifted  to  a  new  vegetation  type.  The  authors 
of  this  model  concluded  that  it  is  the  combination  of  grazing,  drought  or  fire  sup¬ 
pression,  and  military  training  that  makes  it  difficult  for  the  vegetation  to  recover 
from  military  disturbances.  Our  goal  was  to  provide  an  alternative  modeling 
approach  using  an  individual-plant-based  model  (ECOTONE)  developed  for  arid 
and  semiarid  ecosystems  (Peters  2002). 
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2  TASK  1:  STUDY  SITES 

The  first  task  was  to  determine  similarities  in  major  plant  communities  and 
soil  properties  between  Ft.  Bliss  and  the  JER  and  to  quantitatively  identify  those 
plant  communities  that  are  most  similar  between  sites. 

Identification  of  site-level  data 

Identifying  site-level  data  was  not  easy  at  either  site,  although  we  had  more 
support  for  finding  data  on  the  Jornada  Experimental  Range.  A  list  of  long-term 
data  sets  from  the  JER  was  compiled  by  D.P.C.  Peters  in  connection  with  other 
projects  (http://usda-ars.nmsu.edu/data-info/data_index.htm).  For  each  site,  we 
were  interested  in  identifying  soils,  vegetation,  and  climate  information. 

Jornada  Experimental  Range 

The  vegetation  types  at  the  JER  are  well  described.  There  is  a  vegetation  map 
that  was  created  in  1998  by  R.P.  Gibbens  (Fig.  2).  It  is  common  in  semi-arid 
areas  that  vegetation  types  are  strongly  dominated  by  just  one  species.  Based  on 
the  first  dominant  species,  five  major  vegetation  types  have  been  identified:  two 
grassland  types  and  three  shrubland  types  (Fig.  3).  The  two  grassland  types  are 
dominated  by  tobosa  ( Pleuraphis  mutica )  for  playa  grasslands  and  black  grama 
(Bouteloua  eriopoda)  for  upland  grasslands.  The  shrubland  types  are  dominated 
by  honey  mesquite  ( Prosopis  glandulosa),  creosotebush  ( Larrea  tridentata),  and 
tarbush  (Flourensia  cernua).  There  are  a  few  shrublands  at  the  JER  that  are 
dominated  by  a  variety  of  other  species,  but  none  of  these  shrub  species  is  domi¬ 
nant  in  a  very  big  area,  so  vegetation  types  dominated  by  these  species  are  not 
recognized  separately.  Amongst  the  mesquite  shrublands,  R.P.  Gibbens  distin¬ 
guished  three  types  (mesquite  shrublands,  mesquite  dunelands,  and  mesquite 
sandhills)  according  to  the  heights  of  the  sandy  hummocks  formed  around  the 
base  of  the  shrubs.  For  this  project,  we  focused  on  mesquite  shrublands  that  have 
not  yet  formed  significant  hummocks  around  the  base  of  the  shrub,  since  the 
hydrological  implications  of  dune  formation  are  not  completely  understood. 

Soils  information  for  the  JER  is  available  at  the  landscape  scale.  The  soil 
conservation  survey  of  Dona  Ana  county  lists  only  ten  major  soil  associations  for 
the  JER  (Bulloch  and  Neher  1980).  Out  of  these  ten,  four  soil  associations  cover 
the  majority  of  the  area  [the  Onite-Pintura  complex  (OR),  the  Wink-Pintura 
complex  (WP),  the  Dona  Ana-Reagan  Association  (DR),  and  the  Simona- 
Harrisburg  association  (SH)].  The  others  are  OP,  ST,  BJ,  NU,  BK,  WH.  Three  of 
these  minor  soil  types  (NU,  OP,  WH)  cover  the  majority  of  the  area  at  the 
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Figure  2.  Vegetation  map  of  the  Jornada  Experimental  Range  and  the  Chihuahuan 
Desert  Rangeland  Research  Center.  (Unpublished  data  by  R.P.  Gibbens,  figure  by 
D.P.C.  Peters  DACA42-01 -P-01 60.) 
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Figure  3.  Major  vegetation  types  of  the  Jornada  Experimental  Range  and 
the  Chihuahuan  Desert  Rangeland  Research  Center. 

Chihuahuan  Desert  Rangeland  Research  Site  operated  by  NMSU.  We  also  found 
a  more  detailed  map  of  soil  types  at  the  JER  created  in  1962.  This  soil  map  lists 
20  soil  types  for  the  JER.  Unfortunately,  the  description  of  these  soil  types  does 
not  contain  details  for  soil  texture  needed  for  parametrization  of  ECOTONE. 
Additional  information  on  soil  texture  is  available  through  our  collaboration  with 
the  Jornada  LTER  program,  where  detailed  soils  studies  were  conducted  along  a 
transect  running  through  a  bajada  (Wierenga  et  al.  1987).  Furthermore,  the  Desert 
Project  (Gile  and  Grossman  1997)  has  assembled  a  large  amount  of  detailed  soils 
information,  and  B.  Bestelmeyer  (JER-ARS)  is  working  on  correlating  some  of 
these  data  with  vegetation  data  to  develop  a  vegetation-soils  map. 

A  long-term  series  of  climate  data  is  available  from  the  USDA  headquarters 
located  in  the  middle  of  the  Jornada  Experimental  Range  (1915-present).  We  are 
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using  80  years  of  these  historical  weather  data  for  model  parametrization  (191 8— 
1997).  Because  we  were  interested  in  the  effect  of  drought  on  vegetation,  we 
compared  precipitation  during  the  winter  (October-May)  and  the  growing  season 
(June-September)  for  each  year  to  the  long-term  average. 

Fort  Bliss 

Data  from  Fort  Bliss  were  assembled  with  the  help  of  J.S.  Fehmi,  and  they 
were  stored  on  a  CD.  This  CD  contains  data  from  previous  experimental  work  as 
well  as  from  work  conducted  for  the  parametrization  and  verification  of  the 
EDYS  model.  It  is  appended  to  the  second  report  for  the  project  between  NMSU 
and  CERL  (Hochstrasser  et  al.,  in  prep.).  However,  the  data  are  of  limited  use 
since  metadata  are  lacking  and  some  data  are  in  raw  format,  i.e.  they  would  need 
to  be  transformed  for  comparison  with  other  sites.  Soils  information  on  Fort  Bliss 
can  be  obtained  from  the  USDA-National  Resource  Conservation  Service  in  El 
Paso,  TX,  which  has  completed  a  soil  survey  of  Fort  Bliss  and  published  it  on 
line  in  2002  (http://soildatamart.nrcs.usda.gov/Survey.aspx?County=NM013). 

Climate  data  were  obtained  from  the  National  Climatic  Data  Center 
(http://www.ncdc.noaa.gov/oa/ncdc.html)  for  weather  stations  located  near  the 
military  reservation  [El  Paso  International  Airport,  TX  (latitude:  3 1 :49°N,  longi¬ 
tude:  106:23°W,  elevation:  1 194.2  m  asl);  Dell  City,  TX  (latitude:  31:53°N, 
longitude:  105:14°W,  elevation:  1 149.1  m  asl);  and  Orogrande,  NM  (latitude: 
32:23°N,  longitude:  106:05°W,  elevation:  1274.7  m  asl)].  Daily  minimum  and 
maximum  temperature  as  well  as  precipitation  (as  needed  for  model  input)  are 
available  for  El  Paso  Airport  for  1948-2002,  for  Dell  City  for  1980-2002,  and 
for  Orogrande  for  1996-2002. 

Description  and  comparison  of  sites 

Jornada  Experimental  Range  (JER) 

The  Jornada  Experimental  Range  is  a  78,266-ha  area  operated  as  a  ranch  and 
research  station  by  the  USDA  since  1912.  The  site  has  been  moderately  grazed 
by  domestic  livestock.  Over  the  past  century,  numerous  management  techniques 
for  minimizing  shrub  encroachment  have  been  tested  at  this  site,  and  some  prac¬ 
tices  are  still  continued  (e.g.  herbicide  treatment  of  shrubs).  Nevertheless,  shrub 
dominance  increased  from  63%  in  1918  to  92%  in  1998  (R.P.  Gibbens,  unpub¬ 
lished  data).  The  JER  lies  in  an  internal  drainage  basin,  the  Jornada  del  Muerto. 
Elevations  range  from  1,000  to  2,000  m  asl,  with  desert  grasslands  and  shrub- 
lands  occurring  between  1,300  and  1,650  m  asl  (Table  1,  unpublished  data  by  Jin 
Yao).  The  vegetation  is  representative  of  the  northern  Chihuahuan  Desert. 
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The  climate  of  the  Jornada  is  semi-arid  with  mainly  summer  monsoon  rain¬ 
fall.  The  mean  annual  precipitation  is  248  mm  (SD  =  87),  and  the  average 
monthly  temperature  ranges  from  3.8°C  in  January  to  26.1°C  in  July.  Periods 
with  no  rainfall  over  40  days  occur  on  average  about  once  or  twice  a  year 
(Hochstrasser  et  al.  2002a).  The  analysis  of  the  precipitation  pattern  for  the  win¬ 
ter  (October-May)  and  the  growing  season  (June-September)  revealed  the  fol¬ 
lowing  patterns  (Fig.  4):  between  1950  and  1956  and  between  1963  and  1968 
both  the  growing  season  and  winter  precipitation  were  lower  than  the  long-term 
average,  which  resulted  in  severe  droughts  at  the  Jornada  during  which  signifi¬ 
cant  mortality  of  plants  was  observed  (Herbel  et  al.  1972).  Between  1984  and 
1987,  precipitation  was  higher  than  the  long-term  average  in  both  seasons.  Dur¬ 
ing  a  number  of  years  the  winter  was  wet  but  the  summer  was  dry  (using  a  15% 
threshold  of  the  long-term  average)  (1923,  1931,  1979,  1980,  1983,  1992)  (Fig. 
4),  whereas  a  number  of  years  also  had  a  wet  summer  following  a  dry  winter 
(using  a  30%  threshold  of  the  long-term  average)  (1938,  1959,  1974,  1982, 
1996). 


Table  1.  Elevational  range  of  the  five  main  vegetation  types  at  the  Jornada 
Experimental  Range  [Unpublished  data  from  Jin  Yao  (08/13/03)]. 


Dominant 

species 

Minimum 

elevation 

Maximum 

elevation 

Mean 

elevation 

Number  of  30- 
x  30-m  cells 

Boutetoua 

eriopoda 

1313.2 

1634.8 

1361.6 

11,274 

Pleura  phis 
mutica 

1308.2 

1492.1 

1330.5 

5,143 

Larrea 

tridentata 

1352.1 

1530.5 

1368.6 

10,247 

Prosopis 

glandulosa 

1319.8 

1613.1 

1349.0 

12,285 

Flourensia 

cernua 

1306.4 

1490.2 

1329.0 

10,785 
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Decade 


a.  Years  with  15%  above/below  average  precipitation. 


1918-  1928-  1938-  1948-  1958-  1968-  1978-  1988- 
27  37  47  57  67  77  87  96 


■  Oct  -  May 

■  June  -  Sep 

□  Oct  -  May 

□  June  -  Sep 


Decade 


b.  Years  with  30%  above/below  average  precipitation. 


Figure  4.  Comparison  of  the  climate  patterns  in  each  season 
(winter:  October-May,  growing  season:  June-September)  with 
the  long-term  average  precipitation.  The  number  of  years  that 
were  above  (blue)  and  below  (red)  the  long-term  average  is 
depicted  for  each  decade. 
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Fort  Bliss 

Fort  Bliss  military  reservation  comprises  an  area  of  445,154  ha  in  northwest¬ 
ern  Texas  and  southern  New  Mexico.  It  has  been  operated  by  the  military  since 
1889.  The  Bureau  of  Land  Management  is  leasing  parts  of  Fort  Bliss  for  grazing. 
Military  training  occurs  on  the  reservation.  Fort  Bliss  includes  many  different 
habitats  because  of  the  size  of  the  installation  and  the  range  of  elevation  (1,1 64- 
2,692  m  asl).  This  diversity  includes  Chihuahuan  Desert  shrublands  and  grass¬ 
lands  of  several  types  and  Rocky  Mountain  shrublands  and  montane  forests 
(http://www.npwrc.usgs.gov/resource/othrdata/chekbird/r2/fortblis.htm).  Seven 
major  vegetation  types  have  been  distinguished:  Pinyon  juniper,  Mesa  grassland, 
Black  grama  grasslands,  Whitethorn  acacia,  Creosotebush,  Sandsage,  and  Mes¬ 
quite  (Pidgeon  2000).  However,  a  GIS  vegetation  map  from  the  installation 
shows  only  a  gradation  from  lowland  to  upland  vegetation  (Fig.  5). 


Legend 

blissveg 

Value 

High  :  34 


Low :  0 


Figure  5.  Vegetation  map  of  Fort  Bliss. 
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Soils  information  on  Fort  Bliss  is  very  detailed  and  up  to  date  because  the 
NRCS  completed  a  new  soil  survey  in  2002  (Fig.  6).  The  habitat  diversity  of  Fort 
Bliss  sustains  a  large  number  of  plants  and  animals.  Over  300  species  of  birds  are 
found  in  the  various  habitats  on  Fort  Bliss  (http://www.npwrc.usgs.gov/resource/ 
othrdata/chekbird/r2/fortblis.htm).  Notable  mammals  include  the  mountain  lion, 
bobcat,  ringtail  cat,  badger,  bat,  deer,  and  pronghorn  antelope. 

The  climate  of  Fort  Bliss  is  similar  to  the  Jornada  Experimental  Range  in  that 
it  receives  predominantly  summer  rains  (Fig.  7).  The  main  difference  between  the 
Fort  Bliss  climate  and  the  Jornada  climate  is  that  minimum  temperatures  are  gen¬ 
erally  lower  at  the  Jornada  (Fig.  8).  Variation  in  climate  exists  between  different 
parts  of  Fort  Bliss. 


Figure  6.  Soils  map  of  Fort  Bliss  (USDA-NRCS  2002). 
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a.  Mean  monthly  maximum  temperature. 


b.  Mean  monthly  minimum  temperature. 


□  El  Paso 
■  Dell  City 

□  Orogrande 

□  Jornada 


Figure  8.  Comparison  of  weather  data  from  Fort  Bliss  and  the  JES. 
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Conclusions 

Because  of  the  lack  of  detailed  information  about  Fort  Bliss,  it  is  difficult  to 
compare  the  two  sites.  Flowever,  the  climate  is  comparable  between  the  two  sites. 
More  vegetation  types  are  found  at  Fort  Bliss  than  the  Jornada,  but  prevalent 
vegetation  types  at  the  Jornada  can  also  be  found  at  Fort  Bliss. 
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3  TASK  2.  ECOLOGICAL  CHARACTERISTICS  OF 
DOMINANT  SPECIES 

The  second  task  was  to  conduct  a  search  of  existing  data,  peer-reviewed  lit¬ 
erature,  and  gray  literature  to  determine  available  information  on  plants  and  soils 
identified  in  Task  1  and  to  assemble  the  results  into  tables  for  each  plant  species. 

Methods 

We  searched  online  electronic  databases  available  at  New  Mexico  State  Uni¬ 
versity  (NMSU)  and  through  the  USDA-ARS  (AGRICOLA,  BIOSIS,  Current 
Contents)  to  compile  a  bibliography  for  the  five  dominant  species  at  the  JER; 
these  species  are  also  important  at  Fort  Bliss.  The  bibliography  was  compiled  in 
ProCite,  a  bibliography  program  that  allows  flexible  formatting  of  references  and 
that  is  commonly  used  by  ecologists.  Furthermore,  the  Jornada  publications 
archive  and  the  Jornada  LTER  bibliographic  database  were  searched  for  refer¬ 
ences  not  found  in  the  above  resources.  A  total  of  approximately  800  references 
concerning  the  dominant  species  were  identified.  We  also  conducted  interviews 
with  a  number  of  researchers  who  have  experience  with  the  ecology  of  these  spe¬ 
cies  (e.g.,  E.  Fredrickson,  R.P.  Gibbens,  L.  Abbott,  J.  Anderson).  Other  sources 
of  previous  attempts  to  summarize  ecological  characteristics  of  species  are  also 
available.  For  example,  we  found  extensive  ecological  information  on  a  number 
of  arid  lands  species  at  http://www.fs.fed.us/database/feis/. 

Because  the  focus  of  this  project  was  on  species  reactions  to  disturbances  in 
arid  lands,  we  used  a  previously  completed  bibliography  (Fehmi  et  al.  2001)  that 
contains  general  information  on  species  reactions  to  disturbances  (but  not  spe¬ 
cific  to  the  species  of  interest  in  this  project).  This  information  had  to  be  adapted 
to  the  species  of  interest  to  be  useable  for  this  project.  Recently,  characteristics  of 
plant  species  that  contribute  to  their  species-specific  reaction  to  disturbance  have 
been  studied  and  identified  (McIntyre  et  al.  1999;  Weiher  et  al.  1999).  Among 
these  characteristics,  we  focused  on  those  that  are  easy  to  measure.  For  some  dif¬ 
ficult  to  measure  characteristics,  we  used  correlated  and  easier-to-measure  spe¬ 
cies  attributes  (Weiher  et  al.  1999). 

For  easy  access  to  the  references,  we  developed  a  hierarchical  classification 
system  of  keywords  to  help  researchers  find  information  in  the  compiled  bibliog¬ 
raphy  (Appendix  A).  This  resource  has  long  been  needed  by  ecologists  in  arid 
lands,  since  past  attempts  to  assemble  and  gain  an  overview  of  ecological  char¬ 
acteristics  of  these  species  have  been  limited  by  a  lack  of  resources. 

For  the  description  of  the  relation  of  dominant  species  with  soil  characteris¬ 
tics,  we  relied  on  preliminary  research  results  from  B.  Bestelmeyer  (JER-ARS), 
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who  is  studying  vegetation-soil  associations  in  southern  New  Mexico.  In  par¬ 
ticular,  he  has  gathered  information  on  vegetation  at  sites  where  the  soil  was 
described  by  the  Desert  Project  Soil  Monograph  (Gile  and  Grossman  1997).  This 
preliminary  information  is  the  most  comprehensive  and  best  information 
currently  available  on  vegetation-soil  associations  in  southern  New  Mexico. 

Differences  between  grasses  and  shrubs 

Much  has  been  written  about  ecological  differences  between  woody  plants 
and  grasses  (Belsky  1990;  Wilson  1998;  Peltzer  and  Kochy  2001).  Often  the 
difference  in  their  ecological  characteristics  is  further  accentuated  by  their  direct 
and  indirect  effects  on  the  biotic  and  abiotic  environment  (Schlesinger  et  al. 

1990;  Wilson  and  Agnew  1992).  One  of  the  earliest  works  proposed  that  shrubs 
have  access  to  deep  soil  water  that  grasses  cannot  access  (Walker  and  Noy-Meir 
1982).  This  hypothesis  of  vertical  niche  separation  between  the  two  lifeforms  has 
been  supported  by  experimental  evidence  (e.g.  Jackson  et  al.  1996;  Sala  et  al. 
1997;  Dodd  et  al.  1998).  In  environments  where  shallow  water  is  common, 
woody  plants  have  roots  in  both  shallow  and  deep  soil  layers,  and  the  niche  sepa¬ 
ration  between  grasses  and  shrubs  is  not  as  clear  (Wilson  1998).  However, 
grasses  and  woody  plants  differ  in  other  ecological  characteristics,  such  as  their 
root:  shoot  ratio,  the  amount  of  perennial  tissue,  their  competitive  effects  on  sub¬ 
dominant  plants,  their  phenology,  and  their  reaction  and  susceptibility  to  distur¬ 
bances  (Belsky  1994;  Wilson  1998;  Peltzer  and  Kochy  2001;  Hochstrasser  and 
Peters,  submitted). 

Because  root  competition  is  central  to  ECOTONE,  we  surveyed  the  literature 
on  how  root  structure  and  function  differ  between  grasses  and  shrubs.  In  particu¬ 
lar,  we  were  interested  in  how  root  distribution  changes  over  time  and  how  this 
change  affects  the  competitive  abilities  of  grasses  and  shrubs.  A  change  in  root 
distribution  over  time  can  be  observed  if  there  is  root  turnover  (i.e.  if  old  roots 
die  and  new  roots  grow  in  different  places).  Changes  in  root  distribution  over 
time  therefore  depend  on  factors  that  influence  root  growth  and  mortality.  Root 
growth  is  a  function  of  carbon  supply,  but  environmental  factors  such  as  tem¬ 
perature  and  soil  water  distribution  can  be  important  (Luo  et  al.  1995).  Root  mor¬ 
tality  is  a  function  of  the  remaining  above-ground  biomass,  temperature,  and  soil 
water  (Luo  et  al.  1995).  Generally,  root  turnover  increases  with  temperature  (with 
a  Q10  relationship)  (Gill  and  Jackson  2000).  The  phenological  pattern  of  root 
growth  also  differs  between  grasses  and  woody  plants:  root  production  is  high  at 
the  beginning  of  a  growing  season  under  aspen  forest  and  declines  over  the 
growing  season,  whereas  for  mixed-grass  prairie,  root  production  increases  until 
the  peak  of  the  growing  season  and  then  declines  (Partel  and  Wilson  2002);  roots 
in  prairies  seem  to  live  about  twice  as  long  as  in  forests.  Root  turnover  is  34%  for 
total  roots  in  shrublands,  53%  (30-60%)  for  fine  roots  in  grasslands,  50-250% 
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for  fine  roots  in  forest  (Gill  and  Jackson  2000),  20%  for  fine  roots  in  mixed-grass 
prairie,  and  49%  for  fine  roots  in  aspen  forest  (Partel  and  Wilson  2002).  In  blue 
grama,  root  turnover  may  range  anywhere  from  86%/year  (according  to  produc¬ 
tion)  to  35%/year  (according  to  mortality)  (Gill  et  al.  2002).  The  longevity  of 
roots  depends  on  their  diameter:  the  likelihood  of  mortality  decreased  approxi¬ 
mately  6%  with  a  0.10-mm  increase  in  root  diameter,  controlling  for  the  effect  of 
depth  in  the  soil  profile  (Gill  et  al.  2002). 

The  competitive  ability  of  grasses  and  shrubs  depends  on  root  characteristics. 
In  particular,  differences  in  root  surface  area  and  rates  of  resource  uptake  are 
important  determinants  of  competitive  ability  (Eissenstat  and  Caldwell  1988). 
Specific  root  length  is  ten  times  higher  for  grasses  than  for  trees  (Jackson  et  al. 
1997),  and  total  root  length  is  twice  as  long  in  mixed-grass  prairie  than  in  aspen 
forest  (Partel  and  Wilson  2002).  Furthermore,  the  morphological  and  physiologi¬ 
cal  plasticity  in  response  to  nutrient  enrichment  is  an  important  characteristic  of 
plants  contributing  to  their  competitive  ability.  Generally,  the  architecture  of 
roots  changes  in  nutrient-rich  patches  (herringbone  structure)  to  increase  root 
density  (Fitter  1994).  When  comparing  different  lifeforms,  it  is  important  to  con¬ 
sider  that  plasticity  in  biomass  allocation  and  plasticity  in  root  architecture  are 
not  necessarily  related  (Fitter  1994).  Species  with  large  rooting  areas  are  less  able 
to  grow  roots  selectively  in  high-nutrient  patches,  suggesting  a  trade-off  between 
the  ability  to  explore  large  soil  volumes  and  the  ability  to  exploit  nutrient-rich 
patches  (Campbell  et  al.  1991).  Spatial  and  temporal  resource  partitioning  also 
influence  competitive  relationships  between  lifeforms. 

Species  characteristics  of  the  dominants 

Tobosa  (Pleuraphis  mutica)  and  tarbush  ( Flourensia  cernua)  occur  on  finer- 
textured  soil  than  the  other  species  (creosotebush,  mesquite  and  black  grama)  at 
the  Jornada  Experimental  Range  (Fig.  9).  However,  the  distinction  in  terms  of 
soil  differences  is  not  so  clear  (Fig.  9).  Below  we  summarize  the  main  ecological 
differences  and  similarities  between  the  two  dominant  grasses,  tobosa 
( Pleuraphis  mutica)  and  black  grama  (Bouteloua  eriopoda),  and  the  two  domi¬ 
nant  shrubs,  creosotebush  ( Larrea  tridentata )  and  mesquite  ( Prosopis  glandu- 
losa).  Tarbush  ( Flourensia  cernua)  was  no  included  in  this  comparison  because 
very  little  is  know  about  the  ecology  of  this  species  (Hochstrasser  et  al.  2002b). 

Comparison  of  black  grama  and  tobosa 

Black  grama  and  tobosa  occur  on  very  different  soils:  black  grama  tends  to 
occur  on  upland  sandy  soils,  while  tobosa  occurs  in  heavy  clay  soils,  often  in 
playas  that  receive  run-in  water  (Fig.  9)  (Neuenschwander  et  al.  1975;  Wondzell 
et  al.  1996).  Black  grama  is  relatively  long  lived  (Wright  and  Van  Dyne  1976; 
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Figure  9.  Soil  properties  by  depth  of  five  vegetation  types  at  the  Jornada  named  after  their 
dominant  species:  Bouteloua  eriopoda  (BOER),  Flourensia  cernua  (FLCE),  Larrea  triden- 
tata  (LATR),  Pleuraphis  mutica  (PLMU),  and  Prosopis  glandulosa  (PRGL).  (Soils  data  from 
(Gile  and  Grossman  1997);  vegetation  data  from  B.  Bestelmeyer  (pers.  comm.). 


Canfield  1957),  whereas  tobosa  is  shorter  lived  (Wright  and  Van  Dyne  1976). 
However,  tobosa  vegetation  can  be  very  persistent  through  time  (Gibbens  and 
Beck  1988).  The  roots  of  both  species  can  spread  over  1  m  in  depth  (Gibbens  and 
Lenz  2001).  Both  grasses  have  vigorous  vegetative  reproduction;  black  grama 
produces  stolons,  whereas  tobosa  produces  rhizomes  (Nelson  1934;  Paulsen  and 
Ares  1962).  Reproduction  from  seed  is  limited  in  both  species  (Neuenschwander 
1975;  Neilson  1986).  Black  grama  is  more  palatable  to  livestock  than  tobosa.  The 
latter  is  only  palatable  during  the  growing  season  (Paulsen  and  Ares  1962). 


Comparison  of  mesquite  and  creosotebush 

Mesquite  and  creosotebush  both  expanded  their  range  on  the  Jornada 
Experimental  Range  over  the  past  150  years  (Buffington  and  Herbel  1965).  They 
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both  occur  on  a  wide  range  of  different  soil  types.*  It  is  therefore  difficult  to  link 
their  distribution  to  differences  in  soil  characteristics.  Creosotebush  is  an  ever¬ 
green  shrub,  whereas  mesquite  loses  its  leaves  over  the  winter  (November- 
March).  Mesquite  is  a  nitrogen  fixer.  Both  species  are  long  lived,  although  esti¬ 
mates  of  mesquite’s  lifespan  [around  200  years  (Archer  et  al.  1988)]  are  about 
half  of  creosotebush’ s  lifespan  [430  years  (Miller  and  Huenneke  2000)]. 
Creosotebush  seeds  do  not  have  any  particular  adaptation  for  seed  dispersal,  but 
can  be  transported  by  water  and  wind  and  in  the  fur  of  animals  (Boyd  and  Brum 
1983;  Grover  and  Musick  1990).  In  contrast,  mesquite  pods  are  commonly  con¬ 
sumed  by  livestock,  and  the  germinability  of  the  seeds  is  increased  by  scarifica¬ 
tion  as  they  pass  through  the  digestive  system  of  animals  (Brown  and  Archer 
1989;  Archer  et  al.  1988).  In  both  species,  establishment  from  seeds  occurs  infre¬ 
quently.  Both  species  can  resprout  after  the  above-ground  biomass  is  removed. 
Both  species  have  a  wide  root  system  that  can  reach  to  several  meters  deep 
(Gibbens  and  Lenz  2001).  Creosote  started  growth  within  three  days  of  a  simu¬ 
lated  rainfall,  whereas  mesquite  did  not  green  up  (BassiriRad  et  al.  1999).  A 
more  detailed  description  of  the  ecological  characteristics  of  each  species  can  be 
found  in  Appendix  B. 


Personal  communication,  B.  Bestelmeyer. 
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4  TASK  3.  CONCEPTUAL  MODEL  OF  DISTURBANCES  IN 
SHRUBLANDS  VERSUS  GRASSLANDS. 

The  third  task  was  to  build  a  conceptual  model  and  flow  chart  of  plant  and 
soil  interactions  and  dynamics  for  each  location  using  the  information  from  Task 
2. 

Grasslands  and  shrublands  differ  by  the  kind  and  amount  of  feedbacks 
between  soil  resources  and  plants  (Schlesinger  et  al.  1990;  Wilson  and  Agnew 
1992).  In  grasslands,  plants  are  short  lived,  and  therefore  long-term  storage 
occurs  mostly  in  soil  resources,  whereas  shrubs  accumulate  resources  in  their 
biomass  (Fig.  10).  In  the  absence  of  disturbance,  grass  morphology  and  function 
minimize  resource  loss  from  the  soil,  whereas  the  presence  of  shrubs  leads  to  a 
redistribution  and  loss  of  soil  resources  (Schlesinger  et  al.  1990).  Disturbance 
that  impacts  plants  directly  (e.g.  trampling  or  off-road  vehicle  disturbance) 
accentuates  the  loss  of  resource  from  plants  but  not  from  the  soil.  Grasslands, 
where  most  of  the  resources  are  stored  in  the  soil,  are  therefore  likely  to  be  less 
impacted  by  these  disturbances  than  shrublands.  However,  if  disturbances  lead  to 
erosional  processes,  they  greatly  affect  grasslands  and  interrupt  the  regeneration 
cycle  of  grasses  (Fig.  10).  Furthermore,  drought  and  grazing  have  a  greater 
impact  on  grasses  than  on  shrubs.  Grazing  selectively  affects  grasses,  and 
drought  often  shortens  their  lifespan  (Herbel  et  al.  1972).  Therefore,  to  preserve 
grasslands  it  is  very  important  to  minimize  erosion  and  disturbances  that  selec¬ 
tively  affect  grasses,  such  as  grazing.  For  further  conceptual  models  on  how 
plants  are  affected  by  disturbances,  see  Hochstrasser  et  al.,  in  prep.). 

Grassland  Shrubland 


Figure  10.  Conceptual  model  of  disturbances  in  grasslands  and  shrublands. 
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5  TASK  4.  SIMULATION  MODEL  USED:  ECOTONE 

The  fourth  task  was  to  develop  and  validate  a  simulation  model  to  predict  the 
most  likely  future  plant  communities  using  existing  plants,  high/low/average 
weather  conditions,  high/low/average/no  military  training,  and  soil  slope  and 
aspect  for  selected  sites  in  Task  1. 

An  existing  model  (ECOTONE)  was  modified  for  this  project  (Peters  2002). 
ECOTONE  is  a  gap-dynamics,  individual-based  model  in  that  the  recruitment, 
growth,  and  mortality  of  each  plant  in  a  given  vegetation  type  are  simulated 
explicitly  based  on  life-history  traits  of  the  species  or  functional  group;  10-20 
species  or  functional  groups  are  typically  simulated  in  the  model.  A  previous  ver¬ 
sion  of  this  model  (STEPPE)  was  used  to  study  the  impact  of  small-scale  distur¬ 
bances  in  the  shortgrass  steppe  (Coffin  and  Lauenroth  1988).  More  recently, 
ECOTONE  has  been  used  to  predict  changes  in  vegetation  that  may  occur  in 
New  Mexico  and  Colorado  with  directional  changes  in  climate  (Peters  2002; 
Peters  and  Herrick  2002). 

Description  of  the  model 

ECOTONE  is  a  spatially  interactive,  individual-plant-based,  gap-dynamics 
simulation  model  that  simulates  the  recruitment,  growth,  and  mortality  of  indi¬ 
vidual  plants  on  a  small  plot  (Peters  2002)  (Fig.  11).  Recruitment  and  mortality 
have  stochastic  elements,  and  growth  is  determined  by  competition  for  below¬ 
ground  resources  (currently  water).  Plot  size  is  determined  by  the  resource  space 
associated  with  a  full-size  individual  of  the  dominant  species,  and  it  typically 
ranges  from  0.5  to  2  m2  in  desert  grasslands  and  shrublands.  ECOTONE  can  be 
run  either  with  a  spatially  explicit  connection  between  plots  (dependent  land¬ 
scape,  ltype  =  2  -  runs. in)  or  with  no  spatially  explicit  connection  between  plots 
(independent  landscape,  ltype  =  1  -  runs. in).  If  ECOTONE  is  run  with  a  spatially 
explicit  connection  between  plots,  the  seed  availability  in  any  given  plot  depends 
on  seed  dispersal  to  that  plot.  Input  parameters  include  plant  life  history  traits  as 
well  as  environmental  variables  (daily  precipitation  and  temperature,  soil  texture, 
and  disturbance  characteristics).  A  daily  time-step  soil-water  model  (SOILWAT) 
(Parton  1978)  has  been  incorporated  into  ECOTONE  to  represent  seasonal  varia¬ 
tion  in  plant  growth  as  a  result  of  daily  temperature  and  water  availability  by 
depth  in  the  soil  profile.  Growth  and  mortality  of  plants  occurs  on  a  yearly  time 
step.  Output  includes  species  composition;  plant  density;  above-  and  below¬ 
ground  biomass  and  production  by  species,  size,  and  age  distributions;  and  mor¬ 
tality  statistics.  In  what  follows  we  present  details  on  the  processes  in  ECOTONE 
as  depicted  in  Figure  11.  SOILWAT  is  described  in  Appendix  C. 
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Figure  11.  Flow  diagram  of  ECOTONE. 


Dispersal 

The  dispersal  function  used  in  ECOTONE  was  developed  in  1989  (Coffin 
and  Lauenroth  1989),  mainly  with  a  focus  on  grass  seed  dispersal.  The  dispersal 
function  calculates  the  probability  that  a  seed  can  reach  a  plot,  given  its  distance 
from  a  plot  with  that  species  occurring  on  it,  according  to  an  exponential  decay 
function.  The  user  input  is  the  maximum  seed  dispersal  distance  (SDDIST  - 
species. in).  When  calculating  this  distance,  the  user  should  take  into  account  the 
release  height  and  falling  velocity  of  the  seeds  and  the  average  wind  speed  at  the 
study  site  (Peters  2002)  After  the  probability  of  seed  dispersal  to  a  plot  is  calcu¬ 
lated,  it  is  stochastically  determined  if  the  seed  reached  the  plot  or  not  (i.e.  seeds 
will  either  be  available  or  not  for  recruitment  in  this  year).  If  ECOTONE  is  run 
with  a  spatially  independent  landscape  (ltype  =  1  -  runs. in),  the  seed  availability 
in  the  plot  is  1  (i.e.  seeds  are  always  available). 
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Establishment 

Establishment  in  ECOTONE  is  determined  stochastically.  Each  year,  a  num¬ 
ber  (between  MINREC  and  MAXREC  -  species. in)  of  one-year-old  plants  can  be 
recruited  from  a  user-defined  number  of  species  (  between  mnrsp  and  mxrsp  - 
runs.in).  The  total  number  of  recruits  on  each  plot  in  each  year  cannot  exceed  a 
user-defined  number  (mxrec  -  runs.in).  Each  species  has  a  probability  of  recruit¬ 
ment  that  depends  on  seed  availability  (see  Dispersal)  and  a  user-defined  estab¬ 
lishment  probability  (SDECOF  -  species. in).  The  probability  of  recruitment  is 
normalized  to  1.0  for  all  species  eligible  for  establishment.  Species  are  stochasti¬ 
cally  chosen  for  establishment  according  to  their  probabilities  of  establishment 
until  the  maximum  number  of  species  is  drawn  (or  the  maximum  number  of 
recruits  has  been  reached,  whichever  comes  first).  Each  plant  recruited  is  given 
an  initial  biomass  (corresponding  to  a  one-year-old  plant),  which  is  determined  as 
a  percentage  of  the  biomass  of  the  full-grown  plant  [SDFRAT  -  species. in,  the 
full-grown  biomass  being  input  by  FULBIOA  (above-ground  biomass)  and 
FULBIOB  (below-ground  biomass)  -  species.in]. 

Competition 

Plant-available  water  in  each  soil  layer  (the  soil  layer  structure  of  ECOTONE 
can  be  defined  in  soils. in)  and  each  month  is  simulated  using  soils  and  weather 
data  in  SOILWAT  (transpiration  by  layer  by  month).  This  information  is  then 
passed  to  ECOTONE,  and  plants  compete  in  ECOTONE  for  this  water.  The  pro¬ 
portion  of  water  from  each  soil  layer  in  each  month  that  is  distributed  to  each 
plant  depends  on  its  root  surface  area  in  that  layer  as  well  as  its  phenological 
activity.  This  proportion  is  computed  based  the  current  root  biomass,  the 
biomass-to-surface  ratio  (SAREAC  for  coarse  roots,  SAREAF  for  fine  roots  - 
species.in),  the  fine-to-coarse  root  biomass  ratio  (as  specified  in  speclyr.in),  and 
the  phenological  activity  of  the  plant.  The  latter  is  a  scaling  factor  between  0  and 
1  that  is  determined  based  on  the  optimal  temperature  for  growth  (TMPCOF2  - 
species.in)  and  the  average  temperature  of  the  month  [equation  10  in  Peters 
(2002)].  The  relative  amount  of  active  root  surface  area  in  a  given  layer  for  each 
plant  determines  the  proportion  of  the  water  the  plant  receives,  i.e.  the  competi¬ 
tion  is  symmetric  to  the  active  surface  area  of  the  roots. 

Growth 

The  plant-available  water  that  each  plant  receives  is  summed  for  the  year  and 
is  the  amount  of  water  the  plant  can  use  to  grow  and  maintain  its  biomass.  The 
amount  of  biomass  that  can  be  grown  given  the  water  captured  is  determined  by  a 
biomass  conversion  factor  (WUECNST  -  species.in).  Before  growth  can  occur, 
the  perennial  biomass  of  the  plant  has  to  be  maintained.  The  amount  of  water 
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required  to  maintain  the  perennial  biomass  of  the  plant  corresponds  to  a  fraction 
(MTNMODF  -  species. in)  of  how  much  water  was  needed  to  grow  the  biomass 
initially.  If  the  water  captured  by  the  plant  is  less  than  the  water  required  to 
maintain  its  perennial  biomass,  then  clonal  plants  (CLONTYP  =  1  —  species. in) 
can  reduce  in  biomass,  and  other  plants  are  tagged  for  slow  growth.  If  there  is 
more  water  available  than  the  plant  needs  to  maintain  its  tissue,  then  the  plant 
grows.  The  new  biomass  is  allocated  to  different  plant  compartments  according 
to  user-defined  ratios:  the  root-to-shoot  ratio  (PRATIO  -  species. in)  determines 
how  much  of  the  new  growth  is  allocated  above  vs.  below  ground,  the  leaf-to- 
stem  ratio  (PARATIO  -  species,  in)  determines  how  much  biomass  is  allocated  to 
annual  vs.  perennial  tissue  above  ground,  and  the  fine-root-to-coarse  root  ratio 
(speclyr.in)  determines  how  much  biomass  is  allocated  to  fine  vs.  coarse  roots. 

Root  distribution 

Before  the  new  growth  below  ground  is  allocated  to  fine  and  coarse  roots,  the 
new  biomass  is  distributed  among  layers  according  to  relative  amount  of  roots  in 
each  layer.  The  relative  amount  of  roots  per  layer  is  calculated  according  to  a 
user-defined,  species-specific  root  distribution  function  at  the  beginning  of  each 
ECOTONE  run.  The  user  inputs  for  this  function  are  the  depth  of  maximum  root 
biomass  [RTDEPTH(l)  -  species. in]  and  the  maximum  rooting  depth 
[RTDEPTH(2)  -  species. in]  (Fig.  12).  The  root  distribution  is  calculated  accord- 
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Figure  12.  Root  distribution  of  shallow-rooted  grasses 
(Sunetal.  1997,  Table  7). 
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ing  to  a  linear  function  from  the  soil  surface  to  the  depth  of  maximum  root  bio¬ 
mass  and  according  to  an  allometric  function  from  the  depth  of  maximum  root 
biomass  [RTDEPTH(l)]  to  the  maximum  rooting  depth  [RTDEPTH(2)]  (Sun  et 
al.  1997)  (Fig.  12).  The  user  also  inputs  the  relative  amount  of  root  biomass  at  the 
depth  of  maximum  root  biomass  [BMAXD(l)  -  species. in]  and  the  relative 
amount  of  root  biomass  at  maximum  rooting  depth  [BMAXD(2)  -  species.in]. 

Mortality 

At  the  end  of  each  year  the  entire  annual  and  a  certain  percentage  of  the  per¬ 
ennial  biomass  senesces.  The  percentage  of  perennial  biomass  that  is  lost  is 
defined  for  each  plant  compartment:  turnover  percentage  for  coarse  roots 
(CRTORATE  -  species.in),  turnover  rate  for  fine  roots  (FRTORATE  - 
species.in),  turnover  percentage  of  above-ground  perennial  biomass  (LTORATE 
-  species.in),  and  all  of  the  annual  above-ground  biomass.  The  above-ground 
biomass  lost  from  living  plants  is  added  to  the  litter. 

Mortality  of  plants  can  have  one  of  four  sources:  slow  growth,  size  below 
minimum  size,  old  age,  and  disturbance  (Peters  2002).  The  mortality  due  to  slow 
growth  occurs  after  the  plant  has  been  growing  less  than  expected  [i.e.  less  than  a 
certain  percentage  of  its  size  (SLOINC  -  species.in)]  for  more  than  two  years.  Tip 
to  an  age  corresponding  to  half  of  the  species’  maximum  lifespan  (AGEMAX  - 
species.in),  plants  are  exempt  from  mortality  due  to  age.  After  this  age  the  prob¬ 
ability  of  mortality  increases  linearly  until  the  maximum  age  is  reached.  If  the 
probability  is  greater  than  zero,  then  it  is  stochastically  determined  if  a  plant  dies 
of  old  age  or  not.  The  above-ground  biomass  of  plants  that  died  is  added  to  the 
litter. 

Disturbance 

There  are  two  kinds  of  disturbances  currently  in  the  model:  natural  and  tram¬ 
pling.  The  natural  disturbances  occur  stochastically  at  a  user-defined  frequency 
and  have  a  user-defined  size  (disturbances. in).  All  plants  are  killed  in  natural 
disturbances,  except  if  they  can  sprout.  The  probability  of  sprouting  after  distur¬ 
bance  is  user  defined  (VEGGROW  -  species.in).  If  the  plant  sprouts  after  distur¬ 
bance,  then  it  regains  50%  of  its  above-ground  biomass  and  reduces  its  below¬ 
ground  biomass  to  75%  of  its  size  before  the  disturbance  occurred.  The  trampling 
disturbance  routine  is  described  in  Appendix  D. 

Feedback  to  SOILWAT 

Each  year  a  number  of  plant  parameters  are  passed  to  SOILWAT  to  calculate 
the  plant-available  water  as  a  function  of  vegetation  characteristics  (Fig.  11). 
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When  plants  have  reached  their  maximum  size  for  the  year,  the  cumulative 
above-ground  biomass  (standardized  for  area)  of  all  plants  on  a  plot  is  passed  to 
SOILWAT.  After  mortality  and  turnover,  the  amount  of  litter  is  recorded,  and  at 
the  end  of  the  year  the  root  distribution  in  the  plot  is  averaged  between  plants  of 
the  plot.  The  latter  determines  where  in  the  soil  the  plants  take  up  water. 

Model  development 

The  development  of  the  model  was  done  as  a  group  effort  involving  the  sen¬ 
ior  research  scientist,  the  post-doc,  and  the  programmer  working  on  the  model,  as 
well  as  experts  in  fields  researching  different  processes  built  into  the  model.  We 
had  regular  meetings  for  this  development. 

Understanding  the  model  structure  and  function  is  difficult  and  time- 
consuming  because  of  the  large  number  of  variables  involved  and  multiple  func¬ 
tions  of  these  variables.  Initially,  the  algorithms  of  ECOTONE  had  to  be  tested 
and  the  readability  of  the  code  had  to  be  improved  in  order  to  make  changes  to 
the  model.  Once  this  work  was  completed,  we  developed  a  disturbance  routine 
for  trampling  and  vehicle  disturbances.  The  initial  testing  of  ECOTONE  revealed 
that  there  were  problems  in  the  growth  routine,  the  feedback  between  the  vegeta¬ 
tion  and  soil  module  of  ECOTONE,  and  the  root  distribution  routine.  All  of  these 
problems  were  fixed  for  this  project. 

ECOTONE  simulates  a  vegetation  patch  (i.e.  a  small  grid  of  square  plots  that 
are  the  size  of  individual  plants)  by  keeping  track  of  all  individual  plants  within 
this  patch.  This  can  involve  handling  and  storing  a  large  amount  of  data,  and  it 
imposes  large  demands  on  computation  time.  We  therefore  tested  how  long  it 
would  take  to  run  ECOTONE  with  different  spatial  resolution  for  soil  water 
availability.  We  found  that  the  soil-water  model  (SOILWAT)  coupled  with 
ECOTONE  was  indeed  the  limiting  factor  in  the  speed  of  the  simulation  model 
because  SOILWAT  runs  on  a  daily  time  step.  In  the  current  version  of 
ECOTONE,  we  run  SOILWAT  only  once  for  the  entire  patch  per  year  and  not 
for  each  individual  plot  in  the  vegetation  patch  simulated.  Running  SOILWAT 
once  per  plot  and  per  year  would  have  improved  the  correspondence  of  water 
availability  with  vegetation  characteristics  of  the  plot,  but  the  costs  in  computa¬ 
tion  time  did  not  warrant  this  change  to  the  model. 

Mortality  function 

We  improved  the  mortality  function  to  make  it  more  compatible  with  the  rest 
of  the  model.  The  function  that  determined  the  vulnerability  of  plants  to  mortality 
as  they  get  older  was  based  on  a  mortality  function  that  included  ALL  sources  of 
mortality  (Shugart  1984).  Since  we  have  four  sources  of  mortality  (age,  slow 
growth,  smaller  than  minimum  size,  and  disturbance)  in  the  model,  we  changed 
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the  function  to  reflect  increasing  vulnerability  to  disease  etc.  with  increasing  age 
in  the  plants.  We  also  rewrote  some  of  the  output  routines  to  make  output  from 
the  model  more  compatible  with  standard  database  fonnats. 

Disturbance  routine 

We  developed  a  trampling  disturbance  routine  for  ECOTONE  (Appendix  D). 
Sensitivity  analyses 

We  checked  the  sensitivity  of  the  model  to  the  amount  of  biomass  present  on 
the  plots,  the  activity  of  the  vegetation,  and  several  other  species  characteristics. 
When  checking  the  portioning  between  transpiration  and  evaporation,  we  found 
that  the  amount  of  transpiration  increases  with  the  amount  of  active  biomass  and 
litter  in  the  model  up  to  some  maximum,  after  which  it  decreases  because  most  of 
the  precipitation  is  intercepted  by  the  plant  canopy  and  never  reaches  the  soil 
(Fig.  13).  Transpiration  is  a  critical  driver  in  model  dynamics. 


Above-ground  Biomass  (g) 


- Interception 

-  Evaporation 
- Transpiration 


Figure  13.  Variation  of  transpiration  with  amount  of  above¬ 
ground  biomass  in  the  model. 


We  also  tested  the  temperature  response  curve  in  the  model.  This  function 
drives  the  phenology  of  the  plants  in  the  model,  which  is  crucial  for  plant  growth. 
We  found  that  the  total  amount  of  transpiration  in  a  year  is  not  sensitive  to  the 
phenology  of  the  plant  (i.e.  its  optimum  temperature)  (Fig.  14).  However,  the 
amount  of  transpiration  in  each  month  changes  according  to  the  optimum  tem¬ 
perature  for  growth. 


T ranspiration  (cm)  Temperature  effect 
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The  sensitivity  of  the  model  to  optimum  temperature  was  compared  to  its 
sensitivity  to  other  species  characteristics.  We  found  that  the  model  was  sensitive 
(in  decreasing  order)  to  the  water-to-biomass  conversion  factor,  optimum  tem¬ 
perature,  fine/coarse  root  proportion,  and  lifespan.  The  model  was  less  sensitive 
to  establishment  probability  and  seed  dispersal  distance  (results  not  shown). 
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Model  parameters 

We  extracted  species  characteristics  from  references  in  our  literature  review 
(Hochstrasser  et  al.  2002b).  We  also  used  existing  reviews  of  species  characteris¬ 
tics  for  information  (http://www.fs.fed.us/database/feis/;  Smith  et  al.  1997; 

Mabry  et  al.  1977;  Simpson  1977).  However,  the  published  information  was 
insufficient  to  fully  parameterize  the  model,  and  we  had  to  rely  on  expert  knowl¬ 
edge  to  complete  the  parameterization. 

Analysis  and  Results 

We  conducted  a  simulation  experiment  using  ECOTONE  to  compare  recov¬ 
ery  time  after  human  disturbance  in  desert  grasslands  and  shrublands,  both  at  the 
individual  plant  level  and  the  population  level. 

Parameterization  of  the  model 

The  model  was  parameterized  for  black  grama  grasslands  and  mesquite 
shrublands  at  the  Jornada  Experimental  Range.  Nine  species  and  subdominant 
functional  groups  were  simulated  for  each  vegetation  type.  The  soil  was  histori¬ 
cally  dominated  by  black  grama  (Bouteloua  eriopoda)  (Gill  and  Grossman  1 997, 
Photo  1).  Recent  shrub  encroachment  has  made  mesquite  ( Prosopis  glandulosa ) 
the  dominant  vegetation.*  We  compared  the  recovery  dynamics  of  black  grama 
grasslands  versus  mesquite  shrublands  for  this  selected  soil. 

Species  parameters  were  derived  from  a  literature  survey  of  the  dominant 
species  (Hochstrasser  et  al.  2002b).  For  subdominant  species,  existing  parame- 
terizations  of  the  model  were  used.  Vegetation  composition  was  matched  with 
field  records  of  the  vegetation  on  the  soil  used  in  this  model.  Plot  size  was  deter¬ 
mined  according  to  the  resource  space  of  a  full-grown  mesquite  (1.0  m2)  and 
black  grama  plant  (0.25  m2). 

Variation  in  the  recovery  time  of  vegetation  following  a  disturbance  can  be 
caused  by  species  differences  relating  to  the  amount  of  initial  damage  (resistance) 
and  the  amount  of  regrowth/recruitment  after  the  disturbance  (resilience). 
Between  lifeforms,  the  position  of  the  perennating  buds,  the  overall  plant  stature, 
and  the  leaf-stem  architecture  are  important  in  explaining  the  resistance  and 
resilience  of  the  vegetation  (Cole  1995).  In  general,  grasses  are  both  more  resis¬ 
tant  and  more  resilient  to  disturbance  than  other  lifeforms  (Yorks  et  al.  1997). 
Following  these  findings,  we  parameterized  the  model  such  that  grass  species 
were  more  resistant  to  disturbance  (i.e.  they  had  higher  maximum  pressure  toler¬ 
ance)  than  other  lifeforms,  shrub  species  were  intermediate  in  resistance,  and 
forbs  were  least  resistant. 


Personal  communication,  B.  Bestelmeyer. 
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The  model  was  driven  by  80  years  of  daily  weather  data  from  the  Jornada 
Experimental  Range  (1918-1997)  in  southern  New  Mexico,  where  black  grama 
grasslands  and  mesquite  shrublands  are  dominant  vegetation  types.  Model  simu¬ 
lations  were  for  200  years,  during  which  the  80-year  sequence  of  historical 
weather  was  repeated. 

Simulation  experiments 

Determining  the  recovery  time  of  vegetation  after  a  disturbance  is  difficult 
because  of  natural  variation  in  the  system.  For  experiments  on  trampling,  it  has 
been  recommended  that  “recovery”  be  defined  as  the  restoration  of  vegetation 
characteristics,  not  to  pre-disturbance  conditions,  but  rather  to  the  current  condi¬ 
tions  of  a  control  plot  (Cole  and  Bayfield  1993).  The  comparison  with  the  control 
plot  takes  into  account  the  variation  in  vegetation  due  to  climate  and  other  factors 
not  associated  with  the  disturbance.  We  used  this  definition  of  recovery  to  evalu¬ 
ate  our  results. 

Individual  plant  level.  To  determine  disturbance  effects  on  individual 
plants,  we  selected  a  dominant  plant  on  a  plot.  We  applied  a  one-time  disturbance 
when  this  individual  had  reached  an  age  corresponding  to  half  its  maximum  life¬ 
span.  We  repeated  this  experiment  with  a  range  of  disturbance  intensities,  from  a 
low-pressure  impact  (e.g.  one  passage  of  foot  traffic)  to  a  high-pressure  impact 
greater  than  the  plant’s  tolerance  (e.g.  one  passage  by  a  heavy  vehicle).  We 
measured  the  time  it  took  the  plant  to  reach  the  biomass  represented  in  the  no¬ 
disturbance  control  run.  Recovery  of  black  grama  was  compared  to  recovery  of 
mesquite. 

Plant  population  level.  We  indirectly  simulated  the  effect  of  disturbance 
size  on  plant  population  recovery  dynamics  by  simulating  the  effects  of  reduced 
seed  availability  on  recruitment  in  large  disturbances.  Vegetation  dynamics  in  a 
vegetation  patch  uniformly  disturbed  at  high  intensity  (provoking  plant  mortal¬ 
ity)  were  analyzed  comparing  unlimited  seed  availability  (small  disturbance) 
with  a  shortage  of  available  seeds  (large  disturbance).  The  simulated  vegetation 
patch  was  5  x  5  m  in  mesquite  shrubland  (25  plots)  and  2.5  x  2.5  m  (25  plots)  in 
black  grama  grassland. 

Results:  Recovery  time  in  grasslands  and  shrublands 

Individual  plant  level.  Significant  differences  were  found  between  the 
recovery  dynamics  of  a  black  grama  grass  and  a  mesquite  shrub  across  a  range  of 
disturbance  intensities  (Fig.  15).  Black  grama  grass  recovered  within  a  year  from 
low-intensity  disturbances  and  grew  bigger  than  the  control  plant  when  the  dis¬ 
turbance  intensity  was  high.  The  latter  result  could  be  explained  by  the  distur¬ 
bance  impact  on  competitors  of  black  grama,  which  allows  the  already  dominant 


Recovery  time  (%)  P  Total  plant  biomass  (g/plant) 
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Total  plant  biomass  over  time  for  black  grama.  b.  Total  plant  biomass  over  time  for  mesquite. 


Pressure  (%  Maximum  Tolerance) 


c.  Recovery  time  for  black  grama. 


d.  Recovery  time  for  mesquite. 


Figure  15.  Model  output:  Recovery  after  a  range  of  disturbance  intensities  of  individual  black 
grama  and  mesquite  plants. 


plant  on  the  plot  to  capture  more  resources  than  were  previously  available.  In 
contrast,  the  recovery  time  of  a  mesquite  shrub  increased  with  the  disturbance 
intensity.  At  high  intensities,  the  shrub  took  up  to  60  years  to  recover  from  the 
disturbance.  Differences  in  the  recovery  of  these  two  species  can  be  explained  by 
the  amount  of  biomass  removed  by  the  disturbance  as  well  as  their  growth  rate. 
Black  grama  loses  less  biomass  during  disturbance  and  has  a  higher  growth  rate 
than  mesquite. 

Plant  population  level.  The  black  grama  population  was  strongly  reduced 
by  the  disturbance  that  killed  all  plants  in  the  simulated  grassland  patch.  After  the 
disturbance,  the  black  grama  population  took  an  average  of  20  years  to  recover 
(Fig.  16a,  b).  In  contrast,  the  mesquite  population  in  the  simulated  shrubland 
patch  took  longer  to  recover  to  the  control  level  (on  average  about  50  years).  In 
the  control  simulation,  mesquite  biomass  started  to  decline  after  year  150  as  the 


Aboveground  biomass  (g/m‘)  Average  aboveground  biomass  (g/nV) 
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Year 

a.  Black  grama,  unlimited  dispersal. 


Year 


- No  limitation - Dispersal  limited  | 

c.  Black  grama,  limited  dispersal 


Year 

b.  Mesquite,  unlimited  dispersal. 


Year 

—  No  limitation  - Dispersal  limited  | 

d.  Mesquite,  limited  dispersal. 


Figure  16.  Recovery  after  a  disturbance  of  high  intensity  (causing  plant  mortality)  at  the  popula¬ 
tion  level.  Unlimited  dispersal  would  occur  at  a  small  disturbance;  limited  dispersal  would  occur 
at  a  large  disturbance. 


population  reached  its  maximum  lifespan  of  200  years.  In  contrast,  the  individu¬ 
als  in  the  disturbed  populations  were  younger  and  maintained  a  high  level  of 
biomass  at  the  end  of  the  simulation  run. 

Dispersal  limitation  after  the  disturbance,  as  may  occur  in  large  disturbance 
patches,  impacted  the  recovery  dynamics  of  both  dominants  (Fig.  16b,  c):  the 
black  grama  population  still  recovered  relatively  rapidly  but  could  not  attain  the 
average  biomass  of  the  plot  without  the  dispersal  limitation.  This  result  indicates 
that  black  grama  populations  depend  on  relatively  high  recruitment  rates  to  maintain 
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their  population  size.  In  contrast,  recovery  in  the  mesquite  population  was  slower 
than  in  a  non-dispersal-limited  situation,  but  the  population  could  still  attain  and 
even  surpass  the  average  biomass  that  it  maintained  without  the  disturbance. 

Discussion 

These  ECOTONE  simulations  allowed  us  to  capture  important  differences  in 
species  characteristics  between  grasses  and  shrubs  and  to  test  how  these  two  life- 
forms  react  to  disturbances  of  varying  intensities  and  size.  Our  results  correspond 
to  previous  experiments  investigating  the  effects  of  trampling  on  different  life- 
forms  (Cole  1995,  Yorks  et  al.  1997,  Herrick  and  Belnap  2002)  but  contradict 
long-term  observations.  Apart  from  the  biological  differences  between  shrubs 
and  grasses  as  simulated  in  this  work,  the  spatial  distribution  of  anthropogenic 
disturbances  may  not  uniformly  impact  all  plants  (as  assumed  in  this  model).  In 
particular,  mesquite  shrubs  are  commonly  avoided  both  by  foot  and  vehicle  traf¬ 
fic,  which  may  explain  why  shrubs  are  usually  found  in  disturbed  situations 
despite  their  low  resistance  and  resilience. 

The  estimates  of  recovery  time  observed  in  this  research  should  only  be  con¬ 
sidered  guidelines  as  to  the  relative  differences  between  grasses  and  shrubs  in 
arid  environments.  Because  long-term  experiments  are  lacking,  it  is  impossible  to 
verify  if  the  modeled  estimates  are  accurate.  Recovery  times  of  vegetation  also 
depend  on  the  timing  of  the  disturbance  in  connection  with  the  phenology  of  the 
plant.  ECOTONE  currently  simulates  vegetation  growth  and  mortality  on  a 
yearly  timestep,  such  that  seasonal  dynamics  could  not  be  taken  into  considera¬ 
tion. 

Disturbances  in  arid  lands  can  also  lead  to  soil  compaction  that  depends  on 
soil  type  and  soil  moisture  content  at  the  time  of  disturbance.  Soil  compaction 
leads  to  a  decrease  in  infiltration  and  a  decrease  in  plant-available  water.  Surface 
compaction  can  also  inhibit  germination  and  establishment  of  plants.  These  fac¬ 
tors  may  alter  vegetation  dynamics  after  disturbance  and  shift  species  composi¬ 
tion.  Additional  processes  not  simulated  in  the  model  may  also  change  recovery 
dynamics.  For  example,  a  reduction  of  vegetative  cover  can  lead  to  wind  erosion. 

Conclusions 

Grasses  recover  more  rapidly  from  disturbance  than  shrubs  because  of  their 
high  resistance  and  resilience,  both  at  the  individual  plant  level  and  the  popula¬ 
tion  level.  Despite  the  high  resistance  and  resilience  of  grass  plants,  they  are 
more  vulnerable  to  disturbance  than  shrubs  because  of  their  dependence  on 
recruitment  for  maintaining  their  population  size.  Additional  processes,  such  as 
the  effects  of  wind  erosion  and  human  behavior,  may  further  enhance  the  advan¬ 
tages  of  shrubs  for  surviving  trampling  and  vehicle  disturbances. 
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6  SUMMARY  AND  RECOMMENDATIONS 

In  this  work,  we  summarized  the  ecological  characteristics  of  dominant 
shrubs  and  grasses  at  the  Jornada  Experimental  Range.  This  information  is  rele¬ 
vant  to  Fort  Bliss  since  these  vegetation  types  can  also  be  found  there.  We  used 
the  fundamental  ecological  differences  between  shrubs  and  grasses  to  test  for  dif¬ 
fering  reactions  to  disturbances.  First,  we  developed  a  conceptual  model  of  these 
differences.  We  also  incorporated  the  insights  from  our  shrub-grass  comparisons 
into  the  simulation  model  ECOTONE  in  order  to  estimate  recovery  times  from 
disturbances  of  different  intensities  and  sizes.  The  results  obtained  from  the 
simulation  model  indicate  that  this  model  is  capable  of  capturing  the  main  differ¬ 
ences  between  grasses  and  shrubs.  However,  future  work  should  focus  on  incor¬ 
porating  landscape-scale  effects  of  disturbances  and  indirect  processes. 

Research  needs  and  recommendations 

The  simulation  model  we  used  for  this  project  has  the  plant- water  relation¬ 
ships  incorporated  into  it,  but  it  lacks  processes  related  to  plant  nutrient  uptake. 
These  processes  may  be  important  in  arid  lands,  especially  during  times  when 
water  is  available.  More  research  on  the  disturbance  effects  on  nutrient  availabil¬ 
ity  is  needed. 

Soil  erosion  by  wind  and  water  is  often  an  important  consequence  of  distur¬ 
bances  in  arid  lands.  The  redistribution  and/or  loss  of  resources  for  plants  due  to 
erosion  can  influence  plant-plant  interactions  and  future  vegetation  dynamics. 
Because  of  the  strong  effects  of  these  processes  on  community  composition,  we 
recommend  incorporating  these  effects  into  ECOTONE. 

The  damage  caused  by  applying  pressure  to  vegetation  is  not  well  under¬ 
stood,  partly  because  it  is  difficult  to  determine  how  much  pressure  has  been 
applied  to  any  given  spot  in  the  vegetation,  partly  because  the  damage  to  the 
vegetation  has  not  been  well  studied.  Even  though  for  agricultural  vehicles  the 
pressure  distribution  under  the  wheels  has  been  extensively  studied,  less  is 
known  about  the  pressure  under  the  tracks  or  wheels  of  military  vehicles.  Fur¬ 
thermore,  it  would  be  helpful  to  know  how  the  pressure  varies  with  the  driving 
speed  and  style  of  the  driver. 

An  alternative  to  learning  more  about  the  damage  resulting  from  military 
vehicles  at  the  plant  scale  may  be  to  integrate  processes  at  a  higher  organizational 
level  (e.g.  a  patch,  a  landscape)  and  search  for  more  simple  relationships  between 
the  amount  of  training  and  damage  to  the  vegetation.  This  will  not  foster  a 
mechanistic  understanding,  but  it  can  help  to  evaluate  the  amount  of  disturbance 
that  can  be  supported  by  a  certain  type  of  vegetation. 
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In  this  work  we  incorporated  a  study  of  the  ecological  characteristics  of  the 
dominant  species.  However,  to  fully  understand  vegetation  dynamics  after  dis¬ 
turbances,  it  is  important  to  incorporate  subdominant  plants  into  the  study.  This 
poses  a  particular  challenge  because  less  is  known  about  the  ecological  charac¬ 
teristics  of  these  plants  than  about  dominants.  Therefore,  further  studies  on  sub¬ 
dominant  plants  and  their  reactions  to  disturbances  are  needed.  This  can  be  par¬ 
ticularly  crucial  when  invasive  species  are  present  in  the  area  and  pose  a  threat  to 
native  communities. 

Management  implications 

Even  though  shrubs  are  vulnerable  to  above-ground  disturbances,  such  as 
those  caused  by  military  training,  they  are  more  resistant  to  erosional  processes 
than  grasses.  Therefore,  it  is  very  important  to  prevent  wind  erosion  whenever 
possible.  This  can  be  done  by  choosing  areas  that  are  less  vulnerable  to  wind  ero¬ 
sion  (because  of  either  their  soil  characteristics  or  their  high  vegetation  cover)  for 
military  training  (Gillette  and  Adams  1983).  In  contrast,  grasses  are  more  resis¬ 
tant  to  disturbances,  but  because  of  their  high  turnover  rates  and  their  dependence 
on  soil  resources,  they  are  difficult  to  restore  once  the  grass  cover  is  lost.  There¬ 
fore,  it  is  best  to  avoid  the  combination  of  grazing  with  military  disturbance. 
When  grass  cover  is  healthy,  grasses  likely  can  outcompete  shrub  seedlings 
(Brown  and  Archer  1989). 

Plants  also  differ  in  their  vulnerability  to  disturbance  according  to  their  size. 
In  particular,  shrub  seedlings  are  more  vulnerable  to  disturbances  than  adult 
shrubs,  which  suggests  that  if  shrubs  get  established,  it  is  important  to  remove  the 
seedlings  as  soon  as  possible,  while  they  are  still  susceptible  to  control  measures. 
Once  established,  shrubs  contribute  to  resource  redistribution  and  loss  and  often 
make  grass  establishment  impossible  (Wilson  and  Agnew  1992).  We  recommend 
any  measures  that  can  keep  the  lifeform  composition  such  that  grasses  are 
prevalent.  The  structural  changes  in  the  ecosystem  with  a  shift  in  the  dominant 
lifeform  from  grasses  to  shrub  affect  many  ecosystem  processes,  including  the 
resistance  of  the  vegetation  to  future  disturbance. 
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APPENDIX  A.  KEYWORDS  FOR  LITERATURE  DATABASE  ON 
DOMINANT  SPECIES  AT  THE  JORNADA  EXPERIMENTAL 
RANGE 

These  keywords  were  used  for  the  ProCite  database  (named  SPECIES) 
assembled  by  Tamara  Hochstrasser,  2002-2003.  Since  the  objective  of  this  data¬ 
base  was  to  gather  ecological  information  on  dominant  plant  species  at  the  Jor¬ 
nada  Experimental  Range,  the  classification  of  keywords  was  done  from  the 
perspective  of  plant  ecology,  i.e.  categories  get  broader  as  the  focus  shifts  away 
from  the  plant.  For  some  keywords  we  entered  alternative  keywords  into  this 
document  to  show  how  references  were  grouped.  When  this  is  the  case,  the  key¬ 
word  used  is  underlined. 

The  purpose  of  the  assignment  of  keywords  was  to  create  broad  classes  of 
publications.  To  achieve  this  we  devised  a  hierarchical  system  of  keywords. 
There  are  three  broad  categories  that  have  variable  keyword  entries.  They  are  to 
give  a  general  description  of  the  study  system  and  the  organisms  studied.  These 
three  categories  are: 

Scientific  name,  or  Mosses,  Lichens,  Soil  crust 

Location:  Country/U.S.  state,  Name  of  desert,  Continent  (if  NOT  North 

America),  Study  site  (e.g.  Jornada,  Sevilleta,  Santa  Rita,  etc.) 

Vegetation  type  (as  precisely  as  possible):  e.g.  grasslands,  steppe,  desert 

grasslands 

Apart  from  these  three  broad  categories  with  variable  entries,  there  are  also 
another  five  broad  categories  that  have  fixed  entries.  These  eight  categories  are 
considered  the  “top  level.”  The  five  categories  with  fixed  keyword  entries  are 
Species  characteristics,  Community  characteristics,  Kind  of  study,  Disturbance, 
and  Land  use.  Within  each  of  these  five  categories,  attempts  have  been  made  to 
arrange  fixed  keywords  in  a  hierarchical  fashion.  Each  reference  is  assigned 
keywords  at  the  highest  hierarchical  level  first  (i.e.  from  the  top  level).  If  there  is 
more  precise  information,  keywords  from  lower  hierarchical  levels  are  added 
until  further  resolution  does  not  improve  information.  The  idea  is  to  give  the  user 
a  way  to  roughly  classify  references,  before  more  detail  is  added.  To  prevent  the 
classification  of  references  into  these  broad  categories  from  being  made  unclear, 
a  reference  in  one  broad  category  may  contain  low-level  keywords  from  some 
other  broad  category.  For  example,  a  reference  that  belongs  in  the  broad  category 
of  “Plant  physiology,”  may  also  contain  the  keyword  “soil  water  content.” 

These  broad  classifications  will  help  people  who  are  looking  for  references 
concerning  some  subject  matter  for  which  we  have  not  used  a  keyword.  In  such  a 
case  the  user  can  choose  the  broad  categories  that  fit  the  references  he/she  is 
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looking  for  and  add  their  own  lower-level  keyword  as  they  find  references  rele¬ 
vant  to  the  subject  matter.  References  that  are  very  narrow  in  focus  will  have  a 
keyword  from  one  of  the  top  levels  and  then  the  lower-level  specific  keyword, 
without  all  the  intermediary  categories.  For  example,  if  a  reference  reports  the 
specific  leaf  area  of  a  plant,  it  will  be  classified  in  “Morphology:  Leaf  area 
index”  and  NOT  in  “Morphology:  Plant  architecture:  Surface:  biomass  ratio: 
Leaf  area  index.” 

The  hierarchical  classification  of  keywords  may  cause  trouble,  because  the 
classification  is  obviously  quite  subjective,  and  other  people  may  come  up  with  a 
different  hierarchy.  When  looking  for  a  specific  keyword,  it  is  therefore  recom¬ 
mended  to  use  this  list  electronically  such  that  keywords  can  be  identified 
regardless  of  their  hierarchical  level.  The  same  is  true  for  the  ProCite  database, 
which  should  be  able  to  search  for  specific  keywords  (as  long  as  they  are  in  the 
list),  regardless  of  their  hierarchical  level. 


I.  SPECIES  CHARACTERISTICS 

1.  Morphology 

2.  Life  History,  Demography 

3.  Physiology,  Growth 

4.  Chemistry 

5.  Genetics 

6.  Effects  on  Environmental  Fluxes 

7.  Distribution 

8.  Mutualism 

1.1.  Morphology 

1.1.1.  Lifeform 

Raunkiaer  lifeform 

(Position  of  dormant  buds) 

1.1.2.  Plant  size 

Biomass 

above-ground  biomass 
below-ground  biomass 
total  biomass 
biomass  of  seedling 

Size  distribution 

within  population 
age-size  distribution 
size-biomass  relationship 

1.1.3.  Plant  architecture 

Canopy  architecture 

branching  angles 
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plant  height 
leaf  morphology 
Root  distribution 

depth  of  maximum  root  biomass 
maximum  rooting  depth 
root  biomass  distribution  function 
Surface:  biomass  ratio 
leaf  area  index 
specific  leaf  area 
root  surface:biomass  ratio 
Size  of  different  plant  organs 
root:  shoot  ratio 
leaf:  stem  ratio 
reproductive :  vegetative 
fme:coarse  root  ratio 

1.1.4.  Seed  morphology 

Seed  mass,  seed  weight 
Seed  shape 
Seed  surface 

1.1.5.  Plasticity 


1.2.  Life  History,  Demography 

1.2.1.  Regeneration 

Reproductive  age 
Pollination 

pollination  mode 
pollination  efficiency 
Dispersal 

dispersal  mode 
dispersal  distance 
Seed  bank 

spatial  distribution 
seed  bank  persistence 
Germination 

minimum  water  potential  for  germination 
temperature  for  germination 
time  of  germination 
Establishment 

(seedling  survival) 
germinatiomestablishment  ratio 
time  of  establishment 
weather  conditions  for  establishment 
Vegetative  reproduction 

importance  of  vegetative  reproduction 
Sprouting 

amount  of  resprouting 
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1.2.2.  Mortality 

(see  also  DISTURBANCE) 

Lifespan 

Plasticity 

Timing  of  mortality 
/.  2. 3.  Population  dynamics 
Reproductive  rate 
Life  stages 
Growth  rate 
Mortality  rate 


1.3.  Physiology,  Growth 

1.3.1.  Phenology 

Leaf  production 
Root  production 

fine  root  production 
coarse  root  production 
Flowering 
Seed  production 
timing 

amount  of  seed  production 
Leaf  turnover 

leaf  longevity 
leaf  age 
timing 
Root  turnover 

total  root  turnover 
fine  root  turnover 
coarse  root  turnover 
Timing  of  seed  release 
I.  3. 2.  Carbon  fixation 

Primary  productivity 
Photosynthesis 

photosynthetic  pathway 
optimal  temperature  for  growth 
effect  of  shading 
photoperiod 
Respiration 
Transpiration 
WUE 

water  uptake 
plant  water  potential 
1.3.3  Nutrient  uptake 

Nitrogen  fixation 
Nitrogen  uptake 
NUE 

Phosphorus  uptake 
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Potassium  uptake 

1.3.4.  Allocation 

Carbon  allocation 

carbon  storage 
carbon  transport 
Nitrogen  allocation 

nitrogen  reabsoiption 
nitrogen  transport 
Plant  hormones 
1.3.5  Trade-off 

1.4.  Chemistry 

1.4.1.  Biomass  composition 

C:N  ratio 

Lignin  content 

Carbon  content 

Nitrogen  content,  amino  acids 

Litter  composition 

1.4.2.  Secondary  compounds 

1.4.3.  Caloric  content 

1.5.  Genetics 

1.5.1.  Ecotypes 

1.5.2.  Population  level  variability 

1.5.3.  Species  level  variability 

1.6.  Effects  on  Environmental  Fluxes 

1.6.1.  Effects  of  canopy 

Interception 

Shading 

1.6.2.  Effects  of  litter 

1.6.3.  Effects  of  roots 

1.7.  Distribution 

1.7.1.  Range 

1.7.2.  Climatic  constraints 

1. 7.3.  Soil  constraints 

Soil  texture 
Caliche 
1.7.4  Elevation 

1.7.5.  Spatial  pattern 

1.7.6.  Aspect 

1.8.  Mutualism 

1.8.1.  Mycorrhizae 

1.8.2.  Plant-animal  interaction 
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II.  COMMUNITY  CHARACTERISTICS 

1.  Vegetation  structure 

2.  Soil  properties 

3.  Animal  community 

4.  Microsites 

5.  Plant-plant  interaction 

6.  Carbon  dynamics 

7.  Water  dynamics 

8.  Nutrient  dynamics 

9.  Light 

10.  Biotic-abiotic  feedbacks 

11.1.  Vegetation  structure 

11.1.1.  Community  composition 

Cover 

Density 

Species  richness 
Functional  groups 
Species  list 
Exotics 

11.1.2.  Succession 

State  and  transition  model 
Shrub  encroachment 
Ecotone 


11.2.  Soil  properties 

11.2.1  Surface  soil 

11.2.2.  Soil  profile 

11.2.3.  Geomorphology 


11.3.  Animal  community 

11.3.1.  Micro-arthropods 

11.3.2.  Rodents 


11.4.  Microsites 

11.4. 1.  Plant  coexistence  pattern 

11.4.2.  Micrometeorology 

Temperature 
Soil  moisture 

11.4.3.  Microtopography 
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11.5.  Plant-plant  interaction 

11.5.1.  Competition 

Water  relations 
Allelopathy 

11.5.2.  Facilitation 


11.6.  Carbon  dynamics 

11.6. 1.  Amount  of  litter 
II.  6. 2.  Litter  decomposition 
Buried 

fungal 

microbial 

Surface 

termites 

abiotic  decomposition 

11.6.3.  Carbon  distribution 

Microbial  biomass 

11.6.4.  Carbon  cycling  in  soil 

Microbial 

Fungal 

11.6.5.  Carbon  dioxide  fluxes 


11.7.  Water  dynamics 

11.7.1  Evapotranspiration 
Potential 
Actual 

11.7.2.  Water  balance 

Actual  evaporation 

(for  Transpiration  see  Species  Characteristics  ->  Physiology) 

Run-off 

(for  Interception  see  Species  Characteristics  Effects  on 
environmental  fluxes) 

Infiltration 

Soil  water  dynamics 

11.7.3.  Soil  water  content,  Soil  moisture  content 

Content 

Soil  water  potential 

Soil  moisture-soil  water  potential 

Groundwater  table 

11.7.4.  Climate  data 

Precipitation 

Temperature 

soil  temperature 

Relative  humidity,  vapor  pressure  deficit 
Wind  speed 
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11.8.  Nutrient  dynamics 

11.8.1.  Nitrogen 

Nitrogen  distribution 
Nitrogen  cycling  in  soil 
Ammonium  volatilization 

11.8.2.  Phosphorus 

Phosphorus  distribution 
Phosphorus  cycling 

11.8.3.  Potassium 

Potassium  distribution 
Potassium  cycling  in  soil 

11.8.4.  Other  chemical  elements  in  soil 


11.9.  Light 

11.9.1.  Albedo 

11.9.2.  Light  absorption 


11.10.  Biotic-abiotic  feedbacks 

II.  10.1.  Desertification 

11.10.2.  Indirect  effects 


III.  KIND  OF  STUDY 

1.  Observation 

2.  Experiment 

3.  Modeling 

4.  Review 

111.1.  Observation 

Field  survey 
Remote  sensing 
Label 
Long-term 


III. 2.  Experiment 

Greenhouse 
Plant  removal 
Amendment 
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Exclosure 

Isotopes 

Rainfall  simulation,  irrigation 
Remediation 


III. 3.  Modelling 

Deterministic 

Stochastic 


III. 4  Review 


IV.  DISTURBANCE 

1.  Generic  disturbance 

2.  Fire 

3.  Herbivory 

4.  Granivory 

5.  Digging 

6.  Extreme  events 

7.  Erosion 

8.  Dust 

9.  Restoration 


IV.1.  Generic  disturbance 

IV.  1.1.  General  disturbance  characteristics 
IV.  1.2.  Disturbance  history 

IV. 1.3.  Disturbance  interaction  (i.e.  interaction  between  disturbances) 


IV.2.  Fire 


IV.3.  Herbivory 

Clipping 

Insects 

Rodents 

kangaroo  rats 
Lagomoiphs 
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IV.4.  Granivory 

Insects 

Rodents 

kangaroo  rats 

Ants 


IV.5.  Digging 

Rodents 

kangaroo  rats 
Lagomorphs 


IV.6.  Extreme  events 

Drought 

Temperature  stress 

Pulses  in  resource  availability 


IV.7.  Erosion 

Wind,  aeolian 
Water 


IV.8.  Dust 


IV.9.  Restoration 


V.  LAND  USE 


1. 

Grazing 

2. 

Military 

3. 

Roads 

4. 

Brush  control 

5. 

Multiple  land  use 

V.l.  Grazing 

V.l  7. 

Cattle 

Nutrition 

diet  needs 
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diet  composition 
palatability 
Toxicity 
Behavior 
Trampling  effects 
Management 
V.  1.2.  Sheep 
V.1.3.  Goats 
V.  1.4.  Plant  response 

Hormonal  responses 
Chemical  responses 


V.2.  Military 

V.2.1.  Tanks 
V.2.2.  Bombs 
V.2.3.  Infantry 


V.3.  Roads 

V.3.1.  Off- road  vehicles 
V.  3. 2.  Hiking  trails 


V.4.  Brush  control 

V.4.1.  Mechanical 
V.4.2.  Herbicide 
V.4.3.  Techniques 


V.5.  Multiple  land  use 

V.5.1.  Land  use  -  disturbance  interaction 
V.  5. 2.  Land  use  history 

V.5.3.  Interactive  effects  of  different  land  management 
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APPENDIX  B.  SPECIES  CHARACTERISTICS  OF  THE  DOMINANT 
SPECIES  AT  THE  JORNADA  EXPERIMENTAL  RANGE 

The  tables  contain  information  on  the  ecological  characteristics  of  the  domi¬ 
nant  species  of  the  Jornada  as  found  in  the  literature. 

List  of  tables: 

Table  Bl.  Key  references 

Table  B2.  Habitat  of  species 

Table  B3.  Life  history  and  morphology 

Table  B4.  Root  distribution 

Table  B5.  Seed  production,  dispersal 

Table  B6.  Recruitment  characteristics 

Table  B7.  Phenology  and  growth 


Table  Bl.  Key  references. 


Species  name 

References 

Boutetoua  eriopoda 
(BOER) 

1.  Simonin,  K.A.  (2000)  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky 
Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

2.  (84) 

3.  (85,  86) 

4.  (73) 

5.  (69) 

Pleuraphis  mutica 
(PLMU) 

1.  Uchytil,  R.J.  (1988)  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky 
Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

2.  (73) 

3.  (70) 

Flourensia  cernua 
(FLCE) 

1.  Korthuis,  S.L.  (1988)  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky 
Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

Prosopis  glandulosa 
(PRGL) 

1.  Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky 
Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

2.  (83)  p.  167-174 

3.  (77) 

4.  (82) 

Larrea  tridentata 
(LATR) 

1.  Marshall,  K.A.  (1995)  Korthuis,  S.L.  (1988)  U.S.  Department  of  Agriculture, 
Forest  Services,  Rocky  Mountain  Research  Station,  FSL. 
http://www.fs.fed.us/database/feis/ 

2.  (83)  p.  93-106 

3.  (76) 

4.  (92) 

5.  (10) 
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Table  B2.  Habitat  of  species. 


Abb  rev. 

Habitat 

Soil  types 

Elevation 

BOER 

Desert  grasslands* 

Well-drained  sandy  and  gravelly  soilsf 

AZ: 914-1524  m 

CA:  800-1,900  m 

CO:  1,340-1,460  m 

PLMU 

Playa  grasslandsf 

Relatively  impervious  clay  soils** 

AZ: 750-1800  m 

NM:  914-1,981  m 

UT: 488-853  m 

FLCE 

Desert  shrublands** 

Clay  -  clay  loam  soils’1 

400-1,981  m  in 
Chihuahuan  desert 

PRGL 

Desert  shrublandstf 

Wide  variety  of  soils  including  loamy  sand,  sandy 
loam,  calcareous  silt  loam,  noncalcareous  silt 
loam,  gravelly  sand  loam,  deep  sandy  loam,  and 
calcareous  clay*** 

Positively  correlated  with  soil  depth  (65) 

<  1,387  m 

LATR 

Desert  shrublands*** 

Calcareous,  sandy,  and  alluvial  soils  with  underlain 
caliche  hardpan*** 

Absent  from  soils  with  strong  argilic  horizon  (59) 

<  1,515  m  in  all 
deserts 

‘Simonin,  K.  A.  (2000).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station, 
FSL.  http://www.fs.fed.us/database/feis/ 


' Uchytil,  R.  J.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station, 
FSL.  http://www.fs.fed.us/database/feis/ 

“Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station, 
FSL.  http://www.fs.fed.us/database/feis/ 

t+Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station, 
FSL.  http://www.fs.fed.us/database/feis/ 

“‘Marshall,  K.  A.  (1995).  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky 
Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 


Table  B3.  Life  history  and  morphology. 


Abbrev. 

Lifespan 

Lifeform 

Plant  size 
distrib. 

Net  primary 
productivity 

Size/biomass 

relationship 

Full-grown  above-ground  biomass 

Root:shoot 

ratio 

Leaf:  stem  ratio 

BOER 

35-40  years* 

Individual  tillers  seem  to  live 
about  2  years  (43) 

Max.  28  years  (98) 

4%  of  a  cohort  lived  for  14 
years  (31) 

Perennial 

grass 

(75) 

(69) 

38-63  g /  m2/ 
growing 
season  (6) 

(86) 

371.98  g/m2  (33) 

262  g/m2  (79) 

19.8-72.8  g/m2  (depending  on 
cover)  (96) 

378  g /  plant  (85) 

1.45  (79) 

0.3  (86) 

PLMU 

Max.  7  years  (98) 

Pleuraphis  rigida  lives  over 

1 00  years  (24) 

Perennial 

grass 

112.2  g/m2 
(yields)  (70) 

311  g/m2  (51) 

136  g/m2  (37) 

207  g/m2  (1) 

204  g/m2  (1980,  26  g/m2  green),  210 
g/m2  (1981,23  g/m2  green)  (2) 

FLCE 

Perennial 

shrub 

(56) 

164  g/m2  (79) 

PRGL 

1 00-200  yearsf 

Perennial 

shrub 

(81) 

(5) 

17.8  kg/plant/yr 
or  0.78 
kg/m2/yr  (81) 

(14) 

(56) 

258  g/m2  (79) 

6.2  kg/shrub  or  155.2  g/m2  (14) 

219  kg  (big  shrub)  (50) 

43-760  kg/shrub  or  1 .9-8.5  kg/m2  of 
canopy  area  (81) 

0.19  (14) 

0.03  (50) 

0.12,  but  leaf  production 
is  31.6%  of  total 
production  (81) 

(current  growth  & 
leaves)/  (branches  and 
wood)  =  0.17  (Jornada) 

(15) 

LATR 

very  long  lived** 

(Chihuahuan  Desert 
ca.  430  years)  (61) 

Over  100  years  (24) 

11,700  (94) 

1,250  (58) 

Perennial 

shrub 

(75) 

(41) 

(59) 

(13) 

(49) 

(7) 

(9) 

(11) 

(28) 

350-385  kg  dry 
wt  /  ha  (10) 

910  kg/ha/yr 
(California)  or 
917  kg/ha/yr 
(Arizona)  (81) 

(62) 

(29) 

(56) 

60-299  kg/ha  (Mojave)  (13) 

252.56  g/m2  (33) 

204g/m2  (79) 

330  g/m2  (California)  389  g/m2 
(Arizona)  (81) 

120-1 ,820  g  /  plant  (29) 

0.64-6.78 
(mean  2.4)  or 
0.2-2. 7 
(mean  0.9) 

(10) 

0.22-2.5  (76) 
0.33  (44) 

45-70%  of  above-ground 
vegetative  production  is 
allocated  to  leaves  (10) 

0.42  (44) 

*  Simonin,  K.  A.  (2000).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 
f  Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

**  Marshall,  K.  A.  (1995).  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 


Ol 

oi 


Simulation  of  Vegetation  Recovery  from  Military  Disturbances  on  Fort  Bliss 


56 


ERDC/CERL  TR-05-39 


Table  B4.  Root  distribution.  We  did  not  find  any  information  on  fine:coarse  root  ratio,  Fine  root 
surface:  biomass  ratio,  Coarse  root  surface:  biomass  ratio,  or  fine  root  turnover  rate.  This  is  why 
we  used  more  general  information  from  the  literature  as  described  in  the  comparison  between 
grasses  and  shrubs  section. 


Abb  rev. 

Depth  of  maximum 
root  biomass 

Maximum  rooting  depth 

Root  biomass  distribution 
function 

BOER 

5-25  cm* 

1.2  m* 

1 .6  m /  depth  to  caliche  layer  (45) 

1 .22  cm  or  more  (73) 

PLMU 

0-90  cm  (73) 

1 .8  mf(73) 

Over  1  m  (with  caliche)  (45) 

(66) 

FLCE 

>  5  m  (45) 

(66) 

PRGL 

30  cm** 

0-33  cm  (50) 

12  m** 

6-12  m  are  common  (78) 

(50,  66) 

LATR 

20-35  cm+t 

80  cmn 
>  5  m  (45) 

(66) 

*Simonin,  K.  A.  (2000).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research 
Station,  FSL.  http://www.fs.fed.us/database/feis/ 

'Uchytil,  R.  J.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research 
Station,  FSL.  http://www.fs.fed.us/database/feis/ 

**Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research 
Station,  FSL.  http://www.fs.fed.us/database/feis/ 

^Marshall,  K.  A.  (1995).  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky 
Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 


Table  B5.  Seed  production  and  dispersal.  Dispersal  distance  is  an  input  parameter  for  ECOTONE  that  needs  to  be  deduced  from  seed 
mass  and  falling  velocity.  There  is  no  direct  information  on  seed  dispersal  distance  for  these  species. 


Abb  rev. 

Main  regeneration 
mode 

Vegetative 

reproduction? 

Reproductive 

age 

Reproductive 

allocation 

Pollination 

mode 

Amount  of  seed 
production 

Timing  of 
seed  release 

Dispersal 

mode 

Seed  mass 

Seed  shape/ 
falling 
velocity 

BOER 

Stolons*  (69) 

Stolons  (30) 

About  10 

stolons/m2/year 

(69) 

8.3 

stolons/m2/year 

(63) 

2  years* 

51-74  days  after 
germination  in 
green  house 
(natural 
environment  2 
years)  (72) 

October  - 
early 

November 

(69) 

PLMU 

Rhizomest 

In  the  field  few 
seedlings  survive 
(70) 

Rhizomes  (30, 
73) 

Few  fertile  seeds, 
because  of  high 
fungal  infection 
rate  (70) 

FLCE 

Seeds/ 

vegetative 

spread** 

Wind,  Water** 

PRGL 

Seeds' ' 

Seeds  only  (90) 

Sprouting  of 
buried  stems 
(48) 

3  years 

insects 

(bees)* 

Can  be  reduced 
by  drought  or 
freezing  (48) 

Late  spring/ 
early  summer 
(AZ,  NM)  (48) 

Animals* 
Domestic 
livestock, 
rodents,  water 
(48) 

Livestock  (25, 
28) 

LATR 

Vegetative 
growth  / 
rarely  seeds*** 

Mostly  seeds  (48) 

Sprouting  of 
buried  stems 
and  exposed 
roots  (48) 

Shrubs  <  40  cm 
high  did  not 
produce  flowers 
(42) 

After  5-1 3  years 
(48) 

13  years  (32) 

Reproductive 
growth/  total 
growth  0.18  - 
0.8  (35) 

Reproductive 
growth/  total 
growth  0.1  - 
0.8  (40) 

Insects  (74% 
Hymenoptera) 
(23);  pollina¬ 
tion  neces¬ 
sary  for  fruit 
production 
(23) 

1700-  1 '900 
seeds/  shrub 
(23) 

Viable  seeds  at 
dispersal 

represent  only  8% 
of  the  total 
reproductive 
potential  (23) 

Fruits/bush  =  205 
+  66(age  of  shrub 
-  20)  (32) 

Spring  and 
early  fall 
(Mojave  and 
Sonoran)early 
fall  (Chihua- 
huan)  (8,  10) 

Within  1 
month  (in 

July)  94%  of 
seed  were 
dispersed 
(after  bloom 
in  May)  (22) 

Wind,  waterf 
sheet  wash, 
strong  winds, 
caching  of 
seeds  by 
rodents  (48) 

3.5 

mg/seed 

(Sonoran), 

2.8 

mg/seed 

(Mojave), 

1 .6  mg/ 
seed 
(Chihua- 
huan)  (95) 

1.89  mg/ 
seed  (32) 

Threshold 
velocity  for 
wind:  3.9 
km/h(36) 

Simonin,  K.  A.  (2000).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

1  Uchytil,  R.  J.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

1  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

1  Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

1  Marshall,  K.  A.  (1995).  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 
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Table  B6.  Recruitment  characteristics. 


Abbrev. 

Minimum 
water  pot.  for 
germination 

Temp,  range 
for 

germination 

Optimum 
temp,  for 
germination 

Sprouting? 

Seed  bank 
persistence 

Viability 
of  seeds 

Establishment 

probability 

Seedling  biomass 

Root :  Shoot 
seedling  biomass 

BOER 

-1.0  Mpa 
(64) 

20-35°C 

(64) 

2%— 41% 
germinated  in 
different  years 
(69) 

No  seeds  were 
found  at  the 
Jornada  1926 
(52) 

6.4  seedlings/ 
m2/yr  (Arizona) 

(31) 

Seedlings  in  7  out 
of  54  years  (68) 

3.5  months  old  30  g  (16) 

Root:  shoot  of  5  months  old 
plants:  1.2  (38) 

5  weeks  after  development 
of  secondary  roots:  0.4-0. 6 
(39) 

3.5  months  old  0.17  (16) 

PLMU 

up  to  53°C* 

good  under 
39°C  (70) 

39°C 

Yes 

1 0%-40% 
(according  to 
medium) 
germinated  of 
the  seed  crop 
from  1926  (52) 

5  weeks  after  development 
of  secondary  roots:  0.4-0. 8 
(39) 

FLCE 

PRGL 

21-38°C 

(80) 

24- 

27°Ct(89) 
29°C  (80) 

Yes,  with 
3-5  years 
lag  if  dry 

Seeds  remained 
viable  for  10 
years  in  the  soil 
(90) 

42%  germinated 
in  the  field 
(Texas)  (26) 

Over  1 7  years 
about  1  new  tree 
established  in 
every  2m2  (47) 

0.8  g  dry  mass/plant  in 
one  growing  season  (93) 

Around  0.5g  in  4  months 
(74) 

4  months  old  20  g  (16) 

4  months  old  1.41  g  (5) 

0.84-2.42  g /  plant  after 

60  days  according  to  soil 
texture  (27) 

0.6  g  abg,  0.2  g  big  (93) 

Root:shoot  of  4  months  old 
plants  0.675  (74) 

14  weeks  old  0.4-0. 6  (60) 

4  months  old  0.3  (16) 

4  months  old  1.07  (up  to  50 
cm  depth  for  roots)  (5) 

0.18-0.75  according  to  soil 
texture  (27) 

LATR 

needs  >  2.5 
cm  rain** 

After  heavy 
rains  (19) 

Nighttime 
temperature 
>  20°C  (10) 

20-33°C  (8) 

23°C'(8) 

24.5°C  (7) 
late 

summer 

(22) 

Yes1, 
between 
33.3- 
100%  of 
plants  re¬ 
sprouted 
after  fire 
(17) 

No  viable 
meri carps  were 
found  10  months 
after  dispersal 
(22) 

If  stored  at  37°C 

for  30  days 

germinability  of 

seeds 

significantly 

decreased 

(Chihuahuan) 

(7) 

41%  of 
developed 
seeds  were 
viable  (23) 

51% 

(Chihuahuan) 
37%  (Sonoran) 
42%  (Mojave) 

(7) 

45.4%  (61) 
15-76%  (8) 

Low,  especially  in 
dense  grass  cover 
and  creosote  areas 
(97) 

66%  of  seedlings 
(6  months  in 
greenhouse)  (7,  8) 

8,  16  seedlings 
established  (335 
m,  8  yrs)  (18) 

7%/ 3%  over  100 
years  (24) 

1.31  g  (3  mo. 
seedling)(10)  (12) 

After  1-3  years  4.5-15 
cm  in  height  (22) 

4  months  old  10  g  (16) 

0.19(0.22-1.08)  (3  mo. 
seedling)(1 0,12) 

0.33-0.52  (95) 

4  months  old  0.3  (16) 

0.5  (Chihuahuan)  (7) 

*  Uchytil,  R.  J.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 
i  Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

**  Marshall,  K.  A.  (1995).  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 
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Table  B7.  Phenology  and  growth.  Turnover  rates  of  above-ground  biomass  and  estimates  of  the  amount  of  photosynthate  needed  for 
the  maintenance  of  perennial  biomass  were  not  found  in  the  literature.  They  have  to  be  deduced  from  water  use  efficiency  (WUE), 
respiration  rates,  and  growth  rates. 


Abb  rev. 

Phen¬ 

ology 

Optimum  temp,  for 
growth 

Flowering 

Photo¬ 

synthetic 

pathway 

Deciduous? 

Leafing  out 

Growth  rate 

WUE 

BOER 

SU* 

33°C (21) 

August-September 
(5-7  weeks  after 
summer  rains)1  (69) 

Eight  weeks  after 
initiation  of  growth  (43) 

July-September  (after 
rainfall)  (46) 

C4 

remains  green  over 
winter  (43,  73) 

Late  June/  early  July 
depending  on  rainfall  (43) 

First  growth  can  be 
observed  in  late  April/May 
if  year  is  wet  (54) 

Mean  growth  rate  of 
crowns  ranged  between  - 
0.6  -  +1 .4  cm/year  (63) 

437. 7g  H20/g 
biomass(38) 

476  437. 7g  H20/g 
biomass  (73) 

PLMU 

su 

35°C  (34) 

C4 

keeps  growing  over 
the  winter  if 
conditions  are  mild** 

4.73  g/m2/day  (66) 

439. 1g  H20/g 
biomass(38) 

FLCE 

SU 

late  fall,  fruits  mature 
Jan-Marchft 

C3 

Slow 

0.2  mm/day  (66) 

PRGL 

SP-SU* 

29°C  soil  temperature 
(60) 

Below  32°C  (57) 

April-May*** 

Follows  the  same 
pattern  as  leaf 
production,  triggered  by 
photoperiod  (91) 

C3 

loose  leaves 
November  - 
December 
leave  senescence  in 
September  (71) 

loses  leaves  when 
temperature  drop 
below  freezing  (55) 

March-April,  second  cohort 
in  August  (71) 

Soil  tempe-ratures  >  18- 
21  °C  (89) 

End  of  March,  occasionally 
after  summer  rains  (91 ) 

Budbreak  mid-March  - 
mid  April  (Texas)  (4) 

0.1  mm/day  (66) 

annual  percentage 
volume  increase  5.7  % 

(3) 

1,431.7  g  H20/g 
biomass(38) 

LATR 

SP-SU1 

15°C.  Can  be  as  high 
as  25°C  after 
acclimation  (88) 

Photosynthesis  is 
virtually  the  same  in 
January  and  July  (87) 

Different  T  according 
to  season  (67) 

after  rain 

July-August  (42) 

spring  and  late  summer 
(Chihuahuan)  (10) 

after  >  20  mm  rain  and 
a  heat  sum  of  about  443 
degree-days  above 

10°C  (20) 

C3  (10, 

12) 

Evergreen 

Stops  growth  if 
temperature  is 
below  0°C  (55) 
Leaves  live  1  -2 
years  (10) 

Is  able  to  ‘green  up’  within 
10-20  days  after  rainfall 
(53) 

Mainly  in  March-May  and 
August  -  October  (40) 

Spring  after  wet  winter  (1 9) 

Slow3 

Optimum  growth  rate:  4.6 
kg  ha~'  day-1  (55) 
Seedlings:  3%/day  (95) 

Total  abg  growth/  leaf 
biomass  0.01-0.18  (35) 

0.1  mm/day  (66) 

1,628.9  g  H20/g 
biomass(38) 

0.3-1. 5  mg  C02/mg 
H20  (10) 

*  Jornada  Species  List,  Jornada  LTER  program,  http://jornada-www.nmsu.edu/ 

f  Simonin,  K.  A.  (2000).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

**  Uchytil,  R.  J.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

;t  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

***  Steinberg,  P.  (2001).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 

;  :  Marshall,  K.  A.  (1995).  Korthuis,  S.  L.  (1988).  U.S.  Department  of  Agriculture,  Forest  Services,  Rocky  Mountain  Research  Station,  FSL.  http://www.fs.fed.us/database/feis/ 
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APPENDIX  C.  DESCRIPTION  OF  SOILWAT 

This  description  is  based  on  code  (stand-alone  version  1996)  and  Parton 
(1978). 

Introduction 

SOILWAT  is  essentially  what  was  called  the  water- flow  submodel  of  the 
ELM  model,  which  was  a  grassland  simulation  model  for  the  shortgrass  steppe  in 
northern  Colorado.  It  is  a  tipping-bucket  soilwater  model  (Fig.  Cl).  The 
vegetation  parameters  are  community-level  input,  including  total  standing  live 
and  dead  biomass/m2,  litter,  biomass-to-leaf-area  conversion  factor,  and  root 
distribution.  The  soil  characteristics  are  soil  texture  and  water  content  at  field 
capacity  and  at  the  wilting  point. 

The  processes  simulated  on  a  daily  time  step  are  interception,  percolation  of 
water  through  the  soil  (by  fast  and  by  slow  drainage),  evaporation  of  intercepted 
water,  bare-soil  evaporation,  and  transpiration.  Actual  evapotranspiration  is 
calculated  based  on  potential  evapotranspiration  (according  to  Penman’s 
equation)  and  water  potential  in  each  soil  layer.  Actual  evapotranspiration  is 
partitioned  between  transpiration  and  bare-soil  evaporation  based  on  active  leaf 
area.  Transpiration  is  a  function  of  active  leaf  area;  it  occurs  from  all  soil  layers 
where  there  are  roots  present.  Bare-soil  evaporation  can  occur  from  the  upper  soil 
layers  as  specified  by  the  user  (i.e.  the  depth  of  bare-soil  evaporation  can  be 
varied). 

SOILWAT  has  been  used  for  a  long  time  for  numerous  simulation  tasks,  in 
particular  a  number  of  studies  on  germination  and  establishment  probabilities  of 
blue  grama  and  black  grama  (Sala  et  al.  1992;  Lauenroth  et  al.  1994;  Minnick 
and  Coffin  1999;  Peters  2000). 

Detailed  description  of  processes 

The  stand-alone  version  from  1996  shows  slight  differences  from  the  original 
publication  of  the  model  in  1978  (Parton  1978).  In  what  follows  the  stand-alone 
version  is  described.  For  literature  citations,  see  Parton  (1978). 

Interception 

Interception  is  a  function  of  precipitation  and  vegetation  cover  (= 
interception  by  the  standing  biomass)  or  the  amount  of  litter  (=  litter 
interception).  The  equations  are  based  on  Corbett  and  Crouse  (1968),  who 
predicted  interception  for  an  annual  grassland. 
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Figure  Cl.  Flow  diagram  of  SOILWAT. 

It  is  calculated  in  the  following  way: 

Interception  by  the  standing  biomass  (live  and  dead): 

wintstcr  =  parl*0.026*ppt  +  0.094*par2 
if  (wintstcr  >  ppt)  then  wintstcr  =  ppt 
where  ppt  =  precipitation  and  melted  snow  of  the  day 
pari:  if  (vegcov  <  8.5)  then  pari  =  0.9  +  0.04*vegcov 
else  pari  =  1.24  +  (vegcov  -  8.5)*0.35 
par2:  if  (vegcov  <  3.0)  then  par2  =  vegcov*0.333 
else  par2  =  1  +  (vegcov  -  3.0)*0.182 

where  vegcov  (month)  =  pctcover  (month)  *  canopyht  (month) 
where  pctcover  (month)  =  percent  cover  =  biomass/convlai/3 
where  biomass  =  standing  live  and  dead  biomass  of  the  month 
convlai  (month)  =  monthly  amount  of  biomass  needed  to  produce  an 
LAI  of  1,  usually  convlai  =  195.5  (for  arid  systems) 
canopyht  (month)  =  12  +  (34/3.14159)  * 
atan(3. 14159*0. 002*(biomass  -  300)) 
based  on  Conant  (1972). 
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Litter  interception: 

wintlit(NDAY)  =  (0.015*pptleft  +  0.0635)*exp(parl) 
where  pari  =  (-1  +  0.45*alogl0(blitter(month)+l))*alog(10) 
andpptleft  =  ppt  -  wintstcr  (i.e.  precipitation  left  from  precipitation  after 
standing  biomass  interception) 
blitter(month)  =  biomass  of  the  litter  for  the  current  month. 


Water  percolation  in  the  soil 

The  water  is  moved  down  through  the  soil  according  to  the  principle  of  a 
tipping  bucket,  i.e.  each  layer  of  the  soil  is  conceptualized  as  a  bucket,  which  can 
hold  only  the  amount  of  water  that  can  be  held  at  field  capacity.  Any  water 
flowing  into  the  layer  in  excess  of  that  amount  is  moved  to  the  next  deeper  layer 
on  a  daily  time  step  (fast  drainage).  The  water  remaining  in  layers  is  then  slowly 
percolating  into  the  next  deeper  layer  (slow  drainage).  Upward  flow  of  the  water 
occurs  only  through  transpiration  and  evaporation.  Each  layer  maintains  a 
minimum  soilwater  content.  If  water  flow  beyond  the  deepest  layer  is  allowed, 
deep  drainage  occurs;  otherwise,  water  is  accumulating  in  the  deepest  layer. 

Slow  drainage  is  a  function  of  the  difference  between  the  current  water  content  of 
the  soil  layer  and  the  water  content  at  field  capacity,  according  to  the  following 
equation: 


Slow  drainage  =  sdrainpar  *  exp((swc(lyr)-swcfc(lyr))*40/width(lyr)) 


where  sdrainpar 
swc(lyr) 
swcfc(lyr) 
width(lyr) 


slow  drainage  coefficient  (usually  0.06) 
actual  soilwater  content  of  the  layer 
soilwater  content  the  layer  can  hold  at  field  capacity 
width  of  the  layer. 


Water  retention  curve 

For  the  calculation  of  evaporation  and  transpiration  soil  water  potential  is 
used.  Based  on  the  soil  water  content  and  the  texture,  this  is  calculated  according 
to  a  fitted  power  function  from  Cosby  et  al.  (1984): 

swpotent  =  (psis(lyr)/(theta/thetas(lyr))**b/1024 
where  theta  =  soilwater  content  (in  percent  volume)  to  be  converted 
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thetas  =  -14.2*pctsand  -  3.7*pctclay  +  50.5  (i.e.  saturated  water 
content) 

psis  =  10**(-1.58*pctsand  -  0.63*pctclay  +  2.17) 
b  =  -0.3*pctsand  +  15.7*pctclay  +  3.1. 

Potential  evapotranspiration  (PET) 

Potential  evapotranspiration  depends  on  monthly  input,  except  for 
temperature,  which  varies  daily.  Latitude  and  albedo  (fraction  of  light  reflected 
by  the  vegetation)  vary  only  between  sites  for  which  the  model  is  used.  It  is 
calculated  using  Penman’s  equation: 

petfunc  =  ((arads*parA  +  0.27*parB)/(arads  +  0.27))/ 10. 
if  (petfunc<0.01)  then  petfunc  =  0.01 

where  parA  =  0.35*(svapor(avgtemp(jday))  -  fhumid)*(l.  +  0.0098  * 
(windsp( month)  *  24. )) 

parB  =  shwave*(l -reflec)*(0.18  +  0.55*clrsky)  -  ftemp*(0.56  - 
0.092  *  sqrt(fhumid))*(0.10  +  0.90*clrsky) 

arads  =  svapor(avgtemp(jday))*3010.21/(kelvin  *  kelvin) 

svapor  =  saturated  vapor  pressure  for  the  average  temperature  of  the 
day  [calculated  using  the  Clausius-Clapeyron  equation  (Hess 
1959)] 

kelvin  =  avgtemp(jday)  +  273  (i.e.  average  temperature  of  the  day  in 
kelvins) 

windsp  =  average  wind  speed  of  the  month  (user  input) 

clrsk  =  fraction  clear  sky  based  on  the  %  cloud  cover  for  the  month 

shwave  =  short-wave  solar  radiation  on  a  clear  day  according  to  Sellers 
(1965)  (Inputs:  monthly  light  transmission,  Julian  day,  the 
latitude  of  the  site) 

fhumid  =  humidity  (for  the  month)  calculated  using  relative  humidity 
of  the  month  specified  by  the  user  and  saturated  vapor 
pressure  for  the  average  temperature  of  the  day  [calculated 
using  the  Clausius-Clapeyron  equation  (Hess  1959)] 

ftemp  =  ((avgtemp(jday)+273.)*0.01)**4*0.201 

reflec  =  fraction  of  light  reflected  by  the  vegetation  (albedo). 
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Evaporation 

First,  intercepted  water  (by  the  canopy  and  litter)  is  evaporated.  If  PET  is 
bigger  than  the  amount  of  water  intercepted,  all  of  it  evaporates.  PET  is  then 
reduced  by  the  amount  of  water  evaporated  from  interception.  This  constitutes 
atmospheric  demand  for  total  transpiration  and  bare-soil  evaporation. 

Second,  bare-soil  evaporation  is  calculated  as  a  function  of  atmospheric 
demand,  the  inverse  of  the  soil  water  potential  in  layers  from  which  it  can  occur, 
the  total  biomass  (  including  litter),  and  the  active  leaf  area.  The  calculation 
involves  determining  the  total  amount  of  water  lost  from  the  soil  through 
evaporation  and  partitioning  this  total  amount  amongst  soil  layers  from  which  it 
is  evaporating: 

Total  bare  soil  evaporation  loss 

bserate  =  0.5  +  (1. 1/3. 14159)*atan(3. 14159*0. 06*(par2  - 
parl))*pet(jday)*shadeef*fbse 

where  if  (pet(jday)  <  0.2)  pari  =  3.0 

if  (pet(jday)  <  0.4)  pari  =  (0.4  -  pet(jday))*(-10)  +  5 
if  (pet(jday)  <  0.6)  pari  =  (0.6  -  pet(jday))*(-15)  +  8 
if  (pet(jday)  >=  0.6)  pari  =  8 
par2  =  1 5  -  avswp 

avswp  =  average  soil  water  potential  over  all  layers 

~  X(i=l,nelyrs)  width(i)*ecoeff(i)*swpotent/sumecoeff 
swpotent  =  soilwater  potential 

ecoeff(i)  =  evaporation  coefficient  for  the  current  soil  layer 

width(i)  =  width  of  the  current  soil  layer 

sumecoeff  =  sum  of  all  the  evaporation  coefficients 

pet(jday)  =  potential  evapotranspiration  of  the  day 

shadeef  =  (1  -  (totagb/999))  (effect  of  shading  by  the  above-ground 
biomass) 

totagb  =  total  above-ground  biomass,  i.e.  the  sum  of  the  standing  live 
and  dead  biomass  plus  the  biomass  of  the  litter,  if 
(totagb>999)  bserate  =  0 

fbse  =  fraction  of  water  loss  from  bare  soil  evaporation 
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fbse  =  exp(-blivelai*l)  +  0,  if  (fbse  >  0.995)  fbse  =  0.995 
blivelai  =  active  leaf  area. 

This  evaporation  loss  is  then  adjusted  so  that  is  does  not  exceed  the 
atmospheric  demand. 


Actual  evaporation  from  each  soil  layer 


evap  =  (swpfrac(lyr)/sumswp)*bserate 


where  swpfrac(lyr) 
ecoeff(lyr)  = 
swpotent  = 
sumswp  = 
bserate  = 


=  ecoeff(lyr)/swpotent 

evaporation  coefficient  for  the  current  soil  layer 
soilwater  potential 
X(i=l,numlyrs)  swpfrac(lyr) 
total  evaporation  loss  as  described  above. 


Transpiration 

Transpiration  depends  on  soil  texture,  the  amount  of  water  in  the  soil,  the 
evaporative  demand  (PET),  and  a  number  of  vegetation  characteristics.  The 
calculation  of  transpiration  is  analogous  to  the  calculation  of  evaporation:  first, 
the  total  transpiration  loss  is  determined  as  a  function  of  PET,  the  inverse  of  soil 
water  potential  in  different  soil  layers,  root  distribution,  and  the  active  leaf  area. 
Second,  the  transpiration  is  partitioned  by  layer.  Both  the  amount  of  roots  in  each 
soil  layer  and  the  active  leaf  area  are  given  as  a  percentage  (0-1). 

Total  transpiration  loss 

bstrate  =  0.5  +  (1. 1/3. 14159)*atan(3. 14159*0. 07*(par2  - 
par  1 ))  *PET  (j  day)  *fbst 

where  if  (PET(jday)  <  0.2)  pari  =  3.0 

if  (PET(jday)  <  0.4)  pari  =  (0.4  -  PET(jday))*(-10)  +  5 
if  (PET(jday)  <  0.6)  pari  =  (0.6  -  PET(jday))*(-15)  +  8 
if  (PET(]day)  >=  0.6)  pari  =  8 
par2  =  28  -  swpavg 
swpavg  =  min(trwtavg(iregion)) 
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trwtavg(iregion)  = 
swpotent  = 
tcoeff(i)  = 
min/max  = 


sumtcoeff(iregion) 
PET(jday)  = 
fbst  = 
fbse 

hi  i  velai  (month)  = 
biomass 
convlai  = 

pctlive  (month)  = 


d(i  =min,max)  tcoeff(i)*swpotent)/sumtcoeff(iregion) 
soilwater  potential 

percent  root  biomass  for  the  current  soil  layer 

minimum  and  maximum  depth  of  the  given  soil  region 
(iregion),  i.e.  (tcoeff(i)  *  swpotent)  is  summed  over  all  soil 
layers  contained  in  a  region 

=  sum  of  the  transpiration  coefficients  in  the  current  region 

potential  evapotranspiration  of  the  day 

fraction  of  water  loss  from  transpiration  =  1  -  fbse 

fraction  of  water  loss  from  bare  soil  evaporation  (see 
“Evaporation”) 

biomass/convlai*  pctlive 

standing  live  and  dead  biomass  of  the  month 

amount  of  biomass  needed  to  produce  an  LAI  of  1 ,  usually 

195.5  (for  arid  systems) 

percent  live  biomass  (monthly). 


Actual  transpiration  from  each  soil  layer 


trans  =  (swpfrac(lyr)/sumswp)*bstrate 


where  swpfrac(lyr) 
tcoeff(lyr)  = 
swpotent  = 
bstrate  = 
sumswp  = 


=  tcoeff(lyr)/swpotent 

percentage  of  roots  in  the  current  soil  layer 

soilwater  potential 

total  transpiration  loss  as  described  above 
I(i=l,numlyrs)  swpfrac(lyr). 
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APPENDIX  D.  TRAMPLING  DISTURBANCE  ROUTINE 
DEVELOPED  FOR  ECOTONE 


Input  file 

#  trampling. in  -  contains  parameters  related  to  trampling  routine 

# 

#  To  use  this  input  file  set  flag  in  runs. in. 

#  Only  one  trampling  event  can  occur  in  any  given  run  of  the  model. 

#  The  disturbance  is  placed  in  the  middle  of  the  defined  landscape 

#  Disturbance  descriptions  of  different  disturbances  are  included  in  this  file 

#  to  ease  the  switch  between  disturbance  type  in  different  runs. 

# 

# 

#  Variable  descriptions: 

# 

#  name  -  (character)  disturbance  name 

#  distyr  -  (integer)  year  of  occurrence  of  disturbance  (if  0: 

#  disturbance  does  not  occur) 

#  dtwi  -  (real)  width  of  all  tracks  (m) ,  make  sure  this  is  possible 

#  given  tackwi  and  n_passes  below 

#  trackwi  -  (real)  width  of  one  track  (m) 

#  pressure  -  (real)  standing  weight  pressure  (kPa  or  kN*m-2) 

#  n^passes  -  (integer)  number  of  passes 

# 


#name 

distyr 

dtwi 

trackwi 

pressure 

n  passes 

trampling 

50 

0.50 

0.5 

11.5 

1 

#  motorcycle 

0 

0.50 

0.5 

78.0 

1 

#  off-road  vehicle 

0 

1.5 

1.5 

19.6 

1 

#  tank 

0 

1.5 

1.5 

128.3 

1 

# 

# 

#  Species  resistance  to  disturbance 

# 

#  name  -  4  letter  abbreviation  as  in  species. in 

#  MAXP  -  (real)  pressure  necessary  to  damage  all  aboveground  tissue  (kPa) 

#  This  parameter  is  a  relative  ranking  by  lifeform.  It  is  determined 

#  in  relation  to  the  amount  of  damage  occurring  under  a  certain 

#  trampling  regime 

# 

#array  of  1  to  nspec 
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bogr 

MAXP 

250. 

boer 

200. 

arlo 

250. 

spcr 

250. 

stco 

250. 

latr 

150. 

gusa 

100. 

wpfb 

100. 

wafb 

100. 

wags 

100. 

Description  of  the  disturbance  routine 

In  a  user-defined  year,  the  disturbance  routine  places  a  track  of  the  desired 
width  across  the  middle  of  the  simulated  landscape  (Fig.  Dla).  More  than  one 
passage  of  a  person  or  a  vehicle  can  occur.  By  defining  the  disturbance  width,  the 
user  can  define  if  additional  passages  (after  the  first  one)  will  occur  on  the  same 
track  as  the  first  passage  or  immediately  adjacent  to  it.  Subsequent  passages  will 
all  be  adjacent  to  each  other  until  the  disturbance  width  is  filled  (Fig.  Dlb).  In  the 
latter  case  the  area  in  the  middle  will  be  disturbed  a  second  time,  until  all 
passages  have  been  placed  in  the  landscape  (Fig.  Die). 

Each  disturbance  agent  will  exert  a  certain  weight,  calculated  according  to  its 
weight  and  the  size  of  its  “footprint.”  This  weight  is  compared  to  the  weight  that 
the  plant  can  withstand,  and  the  biomass  is  reduced  according  to  the  following 
function  (Fig.  D2): 

%  biomass  reduction  above  ground  =  Pressure  exerted  /  Maxpressure 


If  Pressure  exerted/Maxpressure  <  0.5,  then  %  biomass  reduction  below 
ground  =  0. 


Otherwise  %  biomass  reduction  below  ground  =  2*Pressure 
exerted/Maxpressure  -  1 . 

This  gives  plants  with  a  higher  root:  shoot  ratio  an  advantage  over  plants  with 
a  lower  root:shoot  ratio,  when  the  disturbance  occurs. 

All  the  plants  on  the  plot  are  disturbed,  but  the  effect  varies  between  species 
according  the  pressure  they  can  withstand. 
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a)  One  pass 


Disturbance  width 


Disturbance  width 


Disturbance  width 


J  Track  width 


2  tracks  in  middle 


Figure  D!.  Distribution  of  tracks  in  landscape  as  simulated  by  trampling 
disturbance  routine. 
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Relative  pressure 


—  Reduction  aboveground  Reduction  belowground 


Figure  D2.  Damage  to  a  plant  in  relation  to  the  pressure 
exerted  on  the  plant. 
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